Statistical mechanics of Fofonoff flows in an oceanic basin 
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We study the minimization of potential enstrophy at fixed circulation and energy in an oceanic 
basin with arbitrary topography. For illustration, we consider a rectangular basin and a linear 
topography h = by which represents either a real bottom topography or the /3-effect appropriate 
to oceanic situations. Our minimum enstrophy principle is motivated by different arguments of 
statistical mechanics reviewed in the article. It leads to steady states of the quasigeostrophic (QG) 
equations characterized by a linear relationship between potential vorticity q and stream function 
ip. For low values of the energy, we recover Fofonoff flows [J. Mar. Res. 13, 254 (1954)] that display 
a strong westward jet. For large values of the energy, we obtain geometry induced phase transitions 
between monopoles and dipoles similar to those found by Chavanis & Sommeria [J. Fluid Mech. 
314, 267 (1996)] in the absence of topography. In the presence of topography, we recover and confirm 
the results obtained by Venaille & Bouchet [Phys. Rev. Lett. 102, 104501 (2009)] using a different 
formalism. In addition, we introduce relaxation equations towards minimum potential enstrophy 
states and perform numerical simulations to illustrate the phase transitions in a rectangular oceanic 
basin with linear topography (or ^-effect). 

PACS numbers: 05.20.-y Classical statistical mechanics - 05. 45. -a Nonlinear dynamics and chaos - 05.90.+m 
Other topics in statistical physics, thermodynamics, and nonlinear dynamical systems - 47.10.-g General 
theory in fluid dynamics - 47.15.ki Inviscid flows with vorticity - 47.20.-k Flow instabilities - 47.32.-y Vortex 
dynamics; rotating fluids 



I. INTRODUCTION 

The dynamics of the oceans is extremely complex due 
to nonlinear coupling across many scales of motion and 
the interplay between mean and fluctuating fields [1|. 
Although the oceans can be considered as "turbulent" 
in a classical sense, their dynamics also involves wave- 
like phenomena and coherent structures (vortices) like 
monopoles, dipoles or modons, tripoles... Furthermore, 
despite the permanent action of forcing (e.g. induced by 
the wind) and dissipation, the oceans present a form of 
global organization. This is revealed in the existence of 
strong jets like the Gulf Stream or the Kuroshio Current 
and in the observation of a large-scale oceanic circulation. 
In order to understand the dynamics of the oceans, one 
possibility is to develop numerical codes with increas- 
ing complexity. However, the results can be affected 
by the method used to parameterize the small scales. 
Furthermore, numerical simulations alone do not explain 
the phenomena observed. Therefore, in order to under- 
stand the physical output of such numerical simulations 
it can be useful to consider in parallel simple mathemat- 
ical models that can be studied in great detail. These 
academic models can serve as a basis to develop general 



methods (e.g. statistical mechanics, kinetic theories) that 
can be relevant in more complicated situations. 

Early models of wind- induced oceanic circulation have 
been developed by Stommel [2| and Munk [3] but they 
are based on linearized equations and on an artificial 
concept of eddy viscosity. Alternatively, in a seminal 
paper, Fofonoff [4( neglects forcing and dissipation and 
studies the case of a steady free circulation in a closed 
ocean. He considers quasi geostrophic (QG) flows on the 
/3-plane, a common starting point for many dynamical 
studies in meteorology and oceanography. He further- 
more assumes that the ocean has reached a steady state 
characterized by a linear relationship between potential 
vorticity q = uj + h and stream function ip (here, h de- 
notes the topography and the ordinary /3-effect corre- 
sponds to h = by). Finally, he considers an asymptotic 
regime of low energy and provides a simple analytical so- 
lution representing a westward jet with a recirculation 
at the boundary. This solution is now called Fofonoff 
flow [5]. Numerical simulations starting from a random 
initial condition, in forced and unforced situations, show 
that the system can spontaneously generate a Fofonoff 
flow characterized by a linear q — ip relationship pKLQj. 
However, such a linear relationship is not expected to be 



general and more complex flows with nonlinear q — ip re- 
lationships can also be observed. We must keep in mind 
that FofonofT flows provide just an academic model of 
oceanic circulation with limited applications. 

On the theoretical side, several researchers have tried 
to justify the relevance of Fofonoff flows. In real oceans, 
the flows are forced by the wind and dissipated at small 
scales. Niiler [11] and Marshall & Nurser [12] have ar- 
gued that forcing and dissipation could equilibrate each 
other in average and determine a quasi stationary state 
(QSS) that is a steady state of the ideal QG equations 
(with forcing and dissipation switched off). In these ap- 
proaches, the q — ip relationship is selected by the proper- 
ties of forcing and dissipation and the conditions to have 
a linear relationship are sought. In the case of unforced 
oceans, the justification of a linear q — ip relationship has 
been first sought in a phenomenological minimum enstro- 
phy principle. Bretherton & Haidvogel [13j argue that 
potential enstrophy decays due to viscous effects (see, 
however, Appendix [A]) while the energy and the circu- 
lation remain approximately conserved (in the limit of 
small viscosity). They propose therefore that the system 
should reach a state that minimizes potential enstrophy 
at fixed energy and circulation. This leads to a linear 
q — ip relationship like for Fofonoff flows. 

A linear q — ip relationship can also be justified from 
statistical mechanics. A statistical theory of 2D turbu- 
lence was first developed by Kraichnan [14[ in spectral 
space. It is based on the truncated 2D Euler equations 
which conserve only energy and enstrophy (quadratic 
constraints). In the presence of a topography, Salmon, 
Holloway & Hendershott [15] show that this approach 
predicts a mean flow characterized by a linear relation- 
ship between the averaged potential vorticity (q) and the 
averaged stream function (ip). Another statistical ap- 
proach has been developed by Miller [16j and Robert & 
Sommeria [17j in real space. This theory takes into ac- 
count all the conservation laws (energy and Casimirs) of 
the 2D Euler equation and predicts various q — ip rela- 
tionships depending on the initial conditions. However, 
in real situations where the system undergoes forcing and 
dissipation, the conservation of all the Casimirs is abusive 
and has been criticized by Ellis, Haven & Turkington jl8| 
and by Chavanis pjj [20j. In a recent paper [21], we have 
proposed to conserve only a few microscopic moments of 
the vorticity among the infinite class of Casimirs. These 
relevant constraints could be selected by the properties 
of forcing and dissipation. For example, if we maximize 
the Miller- Robert- Sommeria (MRS) entropy at fixed en- 
ergy E, circulation T and microscopic potential enstro- 
phy Ig ', we get a mean flow characterized by a linear 
q — tp relationship leading to Fofonoff flows. This statisti- 
cal approach also predicts Gaussian fluctuations around 
this mean flow. Furthermore, we have shown that the 
maximization of MRS entropy at fixed energy, circula- 
tion and microscopic potential enstrophy T^' 9 ' is equiv- 
alent to the minimization of macroscopic potential en- 
strophy Y^' 9 ' at fixed energy and circulation. This jus- 



tifies an inviscid minimum potential enstrophy principle, 
and Fofonoff flows, when only the microscopic enstro- 
phy (quadratic) is conserved among the infinite class of 
Casimirs. 

The asymptotic limit treated by Fofonoff [4] corre- 
sponds to small energies E, a limit relevant to oceanic 
situations. In the statistical theory, this corresponds to 
a regime of large positive inverse temperatures /? ^> 1. 
In that case, the Fofonoff solution is the unique (global) 
entropy maximum at fixed circulation and energy. On 
the other hand, Chavanis & Sommeria [22| studied the 
case of a linear q — ip relationship in a rectangular domain 
without topography. For large energies E, corresponding 
to sufficiently negative /?, they report the existence of 
multiple solutions. This leads to interesting phase tran- 
sitions between monopoles and dipoles depending on the 
value of a single control parameter Y/y/E and on the ge- 
ometry of the domain (for example, the aspect ratio r 
of a rectangular domain). For T = 0, the maximum en- 
tropy state is a monopole if r < r c = 1.12 and a dipole 
if r > r c . For r / 0, when r < r c the maximum en- 
tropy state is always a monopole and when r > r Cl the 
maximum entropy state is a dipole for small values of 
Y 2 /E and a monopole for large values of Y 2 /E. The 
approach of Chavanis & Sommeria [22| h as been com- 
pleted recently by Venaille & Bouchet [23|, [24| who used 
a different theoretical formalism and provided a detailed 
discussion of phase transitions and ensemble inequiva- 
lence in QG flows with and without topography. For 
low values of the energy, they recover Fofonoff flows Q 
and for large values of the energy, they obtain geom- 
etry induced phase transitions between monopoles and 
dipoles similar to those found by Chavanis & Sommeria 
[22| . They also emphasize the notions of bicritical point 
and azeotropy. Their approach is well-suited for statisti- 
cal mechanics but it may appear a bit abstract to fluid 
mechanicians. By contrast, the approach of Chavanis & 
Sommeria [22J is simpler. In the present paper, we shall 
extend the formalism of |22[ to the case of flows with a 
topography and study how the series of equilibria is mod- 
ified in this more general context. We first consider an 
antisymmetric linear topography h = by in a rectangu- 
lar domain (like in Fofonoff 's classical study) and then 
generalize the results to the case of an arbitrary topogra- 
phy in an arbitrary domain. We recover and confirm the 
main results of Venaille & Bouchet [23|, [24[ and illustrate 
them with explicit calculations and with synthetic phase 
diagrams. We compute the full series of equilibria (con- 
taining all the critical points of entropy at fixed energy 
and circulations) while Venaille & Bouchet [23|, [24[ focus 
on global entropy maxima. The full series of equilibria 
is useful to show the relative position of the stable and 
unstable branches. Furthermore, for systems with long- 
range interactions, metastable states (local entropy max- 
ima), and even saddle points of entropy, can be long-lived 
and therefore relevant in the dynamics [21] . It is there- 
fore important to take them into account. We also give a 
special attention to the existence of a second order phase 



transition that exists only for a particular value of the 
circulation r* (with r* = for an antisymmetric topog- 
raphy) and study how this phase transition takes birth 
as T —> T* through the formation of a "spike" . Finally, 
we introduce simple relaxation equations that converge 
towards the minimum potential enstrophy state at fixed 
circulation and energy and solve these equations numer- 
ically to illustrate the phase transitions in a rectangular 
oceanic basin with linear topography (or /3-effect). Nu- 
merical integration of these equations with a nonlinear 
topography have been previously performed in [25j and 
the present paper develops the theory required to inter- 
pret the results. 



II. THE QUASIGEOSTROPHIC EQUATIONS 

A. A maximization problem 

We consider a 2D incompressible flow over a topogra- 
phy described by the quasigeostrophic (QG) equations 



Let us consider the maximization problem 
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-Aip + ^ + h, (1) 



where q = uo + tp/R 2 + h is the potential vorticity, h 
the topography, R the Rossby radius, ujz = V x u the 
vorticity, xp the stream function and u = — z x V^ the 
velocity field (z is a unit vector normal to the flow). For 
illustration and explicit calculations, we shall consider a 
linear topography of the form h = by. This term can 
equivalently be interpreted as a "/3-effect" in the oceans 
due to the Earth's sphericity. 

The QG equations admit an infinite number of steady 
states of the form 



/wo> 



(2) 



where / is an arbitrary function. They are obtained by 
solving the differential equation 



A^ + — ^ + /i = /(V), 



(3) 



with ip = on the domain boundary. The QG equations 
conserve the energy 



E 



(q — h) ip dr 
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dr, (4) 



and an infinite number of integral constraints that are 
the Casimirs 



/ 



g(q) dr, 



(5) 



where g is an arbitrary function. In particular, all the 
moments of the potential vorticity T n = J q n dr are con- 
served. The first moment T = J q dr is the potential 
circulation and the second moment T2 = J q 2 dr is the 
potential enstrophy. 



max{5'[g] I E[q] 
q 



e, T[q] = r}, 



(6) 



where E and V are the energy and the circulation and S 
is a functional of the form 



/ 



C(q)dr, 



(7) 



where C is an arbitrary convex function (i.e. C" > 0). 
The critical points of S at fixed E and Y are given by the 
variational principle SS — {38E — aST = where j3 and a 
are Lagrange multipliers. This gives C'{q) = — flip — a. 
Since C is convex, we can inverse this relation to obtain 
q = Ftfil) + a) = /(VO where F(x) = (C / )" 1 (-x). We 
note that q f (ip) = —j3/C"(q). Therefore, a critical point 
of S at fixed E and V determines a steady state of the 
QG equations with a monotonic q — ip relationship that 
is increasing for /3 < and decreasing for /3 > . On the 
other hand, this state is a (local) maximum of S at fixed 
E and T iff 



C"(q)(5q) 2 dr--/3 SqSip dr < 0, 



(8) 



for all perturbations 5q that conserve energy and circu- 
lation at first order. 



B. Its different interpretations 

The maximization problem ((6]) can be given several 
interpretations (see [26| for a more detailed discussion): 

(i) It determines a steady state of the QG equations 
that is nonlinearly dynamically stable according to the 
stability criterion of Ellis et al. [18|. In that case, S will 
be refered to as a "pseudo entropy" [26| . This criterion is 
more refined than the well-known Arnol'd theorems [27J 
that provide only sufficient conditions of nonlinear dy- 
namical stability. We note, however, that the criterion 
([6]) provides itself just a sufficient condition of nonlin- 
ear dynamical stability. An even more refined criterion 
of nonlinear dynamical stability is given by the Kelvin- 
Arnol'd principle [28Ll29j. The connections between these 
different criteria of dynamical stability are reviewed in 
|26j . For the particular choice S = —(1/2) J q 2 dr (neg- 
enstrophy), the maximization problem (|6]) determines 
steady states of the QG equations with a linear q — ip 
relationship that are nonlinearly dynamically stable. 

(ii) The maximization problem (|6]) can be viewed as a 
phenomenological selective decay principle (for —S) due 
to viscosity [2Q, [30| . In the presence of a small viscos- 
ity v — >• 0, the fragile integrals —S decay (see, however, 
Appendix [A]) while the robust integrals E and F remain 
approximately conserved. This suggests that the system 
will reach a steady state that is a maximum of a certain 
functional S at fixed E and I\ If we assume that this 
functional is S = — (1/2) j q 2 dr, we recover the ordinary 



minimum enstrophy principle introduced by Bretherton 
& Haidvogel [13j (however, these authors mention in their 
Appendix that other functionals —S of the form (0, that 
they call generalized enstrophies, could be minimized as 
well). We can also justify an inviscid selective decay 
principle due to coarse-graining [22|, [26j. For an ideal 
evolution (no viscosity), the integrals S[q] = — J C(q) dr 
of the fine-grained PV are conserved by the QG equa- 
tions (they are particular Casimirs) while the integrals 
S[q] = — f C(q) dr of the coarse-grained PV increase (see 
Appendix A of |31j). T hese functionals are called gener- 
alized H- functions [3lLl32J|. This suggests that the system 
will reach a steady state that is a maximum of a certain 
H- function S[q] at fixed E and V. If we assume that 
this functional is S = — (1/2) fq 2 dr, we justify an invis- 
cid minimum enstrophy principle due to coarse-graining 
j22j. However, other generalized H functions could be 
maximized as well [26j. It is important to emphasize 
that these principles are purely phenomenological and 
that they are not based on rigorous arguments. As such, 
they are not always true [26|, [33j . 

(hi) The maximization problem (J6j) provides a neces- 
sary and sufficient condition of thermodynamical stabil- 
ity in the Ellis- Haven- Turkingt on (EHT) approach [18j 
where the Casimir constraints (fragile) are treated canon- 
ically so that they are replaced by the specification of 
a prior vorticity distribution x( cr ) encoding small-scale 
turbulence. Indeed, a vorticity distribution p(r, <j) is a 
maximum of relative entropy S x [p] at fixed circulation 
and energy (EHT thermodynamical stability) iff the cor- 
responding coarse-grained PV field q(r) is a maximum of 
a "generalized entropy" S[q] at fixed E and T [18l I20I. l25j . 
In that case, the generalized entropy S[q] is determined 
by the prior x( cr )- For a Gaussian prior, we find that the 
generalized entropy S\q] = —(1/2) fq 2 dr is proportional 
to minus the macroscopic coarse-grained enstrophy, jus- 
tifying a minimum potential enstrophy principle in that 
context. 

(iv) Since the solution of a variational problem is al- 
ways solution of a more constrained dual variational 
problem (but not the converse) [34[, the maximization 
problem (|6)) provides a sufficient condition of thermody- 
namical stability in the Miller- Robert- Sommeria (MRS) 
approach [16|, LL7j where all the Casimirs are conserved. 
Indeed, a vorticity distribution p(r, a) is a maximum 
of MRS entropy S[p] at fixed energy, circulation and 
Casimirs (MRS thermodynamical stability) if it is a max- 
imum of S x [p] at fixed energy and circulation. However, 
the converse is wrong since the Casimir constraints have 
been treated canonically. According to (iii), we conclude 
that a vorticity distribution p(r, a) is a maximum of MRS 
entropy S[p] at fixed energy, circulation and Casimirs if 
the corresponding coarse-grained PV field q(r) is a maxi- 
mum of S[q] at fixed E and T (but not the converse) [35|. 
In that case, the generalized entropy S[q] is determined 
by the initial condition [26J |. For initial conditions that 
lead to a Gaussian PV distribution at statistical equilib- 
rium, we find that S[q] = —(1/2) j q 2 dr. Therefore a 



minimum of coarse-grained potential enstrophy at fixed 
circulation and energy is a MRS thermodynamical equi- 
librium but the converse is wrong in case of ensemble 
inequivalence. In the MRS approach, a minimum enstro- 
phy principle can be justified only if the microcanonical 
and grand microcanonical ensembles are equivalent |26j . 
(v) The maximization problem (J6]) provides a suffi- 
cient condition of thermodynamical stability in the Naso- 
Chavanis-Dubrulle (NCD) approach [21[ where only a 
few Casimirs are conserved (the reason is the same as 
in (iv)). In that case, the generalized entropy S[q] is de- 
termined by the set of conserved Casimirs. For example, 
if we only conserve the microscopic potential enstrophy 
T f 2 g - = JVdr, we find that S[q] = -(1/2) fq 2 dr. In 
that case, the generalized entropy is proportional to mi- 
nus the coarse-grained enstrophy T^' 9 ' = — fq 2 dr. Fur- 
thermore, in this specific case where the constraints are 
linear or quadratic (energy-enstrophy-circulation statis- 
tical mechanics), it can be proven that the maximization 
of MRS entropy at fixed energy, circulation and micro- 
scopic potential enstrophy T^' 9 ' (NCD thermodynamical 
stability) is equivalent to the minimization of macroscopic 
potential enstrophy T^' 9 ' at fixed energy and circulation 
[21] . This justifies an inviscid minimum potential enstro- 
phy principle, and Fofonoff flows, when only the micro- 
scopic enstrophy is conserved among the infinite class of 
Casimirs. 



C. Relaxation equations 

Some relaxation equations associated with the maxi- 
mization problem (J6| have been introduced in [26| . They 
can serve as numerical algorithms to solve this maximiza- 
tion problem. They also provide non trivial dynamical 
systems whose study is interesting in its own right. 

(i) The first type of equations is of the form 



^ + u.V, = V 



f3(t) 



D[Vq 



C"(q) 



J DVq ■ Vip dr 



Vip 



C"(q) 



dr 



(9) 



(10) 



where D(r,t) > is the diffusion coefficient. The bound- 
ary conditions are Jn = where J = — D(Vq+ cm \ VVO 
is the current and n is a unit vector normal to the bound- 
ary. With these boundary conditions, the circulation is 
clearly conserved. On the other hand, the inverse "tem- 
perature" f3(t) evolves in time according to Eq. ([TO]) so 
as to conserve energy {E = 0). Easy calculations lead to 
the i^-theorem: 

2 



5 = /DC»«„(v s+ ^V* 



dr > 0. 



(11) 



C"(q) 

Therefore, the relaxation equations ([9])- (fTQ|) relax to- 
wards a (local) maximum of S at fixed E and T (see 
[26] for a more precise statement). 



Example: If we take S to be the opposite of the poten- 
tial enstrophy S = —(1/2) J q 2 dr, we get 



dq 

di 



+ u-Vq = V-[D(Vq + /3(i)VV>)] , (12) 



m 



j DVq ■ Vip dr 
' jD(Vtp) 2 dr ■ 



(13) 



This equation monotonically dissipates the potential en- 
strophy (f 2 = -25 = -2 / D(Vq + /?(£) V?A) 2 dr < 0) at 
fixed energy and circulation (T = E = 0) until the min- 
imum potential enstrophy state is reached. If we take 
D constant and R = oo (for simplicity), the foregoing 
equations reduce to 



dq 



m 



u-Vq = D [Aq - 0(t) (q - ft)] , 



S(t) + ±(qh) -T 2 {t) + (qh) 



E 



2E 



(14) 



(15) 



where we have used an integration by parts in Eq. (fT3|) 
to obtain Eq. (fl5j) . In particular, for ft = 0, we have 
f3(t) = S(t)/E == -T 2 {t)/2E. As shown in Appendix 
\B\ these relaxation equations are compatible with the 
"Neptune effect" discovered by Holloway [36| and playing 
an important role in oceanic modeling. 

(ii) The second type of relaxation equations is of the 
form 



dq ,_, 



-D(C'(q) + l3(t)^ + a(t)), 

where (3(t) and a(t) evolve in time according to 

(C'(q)4>)+p(t)^ 2 )+a(t)^)=0, 



(16) 



(17) 



with 



m 



a(t) 



T(il>) - A(2E + (hip)) 



T{^ 2 )-^){2E+(hj>)) 



(21) 



(22) 



These equations monotonically dissipate the potential en- 
strophy (f 2 = -25 = -2D J((7+/3^+a) 2 dr < 0) at fixed 
energy and circulation (E = T = 0) until the minimum 
potential enstrophy state is reached. 



III. THERMODYNAMICS OF FOFONOFF 
FLOWS 

A. The maximization problem 

We shall study the maximization problem 



with 



max{S[q]\E[q]=E, T[q]=r}, 
q 



-\\<?dr 



(23) 



(24) 



E 



l -J(q-h)^dv= l -J{V^) 2 dv, (25) 



/ 



qdr, 



(26) 



(C"(q)}+f3(t)W+o : (t)A = 0, 



(18) 



in order to satisfy the conservation of energy and circu- 
lation (A is the domain area, (X) = J X dr and we have 
assumed D constant for simplicity). We shall consider 
boundary conditions of the form C'(q) + a(i) = on the 
boundary so as to be consistent with the steady state for 
which C'(q) + f3ip + a = in the whole domain (recall 
that i(j = on the boundary). Easy calculations lead to 
the H theorem: 

5 = D f(C'(q) + f3(t)iP + a(t)) 2 dr > 0. (19) 



Therefore, the relaxation equations (fT6|) - (fT8]) relax to- 
wards a (local) maximum of 5 at fixed E and T (see [26j 
for a more precise statement). 

Example: If we take 5 to be the opposite of the poten- 
tial enstrophy 5 = — (1/2) j q 2 dr, the relaxation equa- 
tions reduce to 



dq 



Vq = -D[q + P(t)tl> + a(t)] 



(20) 



q = -A^ + ft. 



(27) 



For simplicity, we assume R — >• +oo but the case of finite 
Rossby radius can be treated similarly and the main re- 
sults are unchanged. As discussed in Sec. Ill Bl the max- 
imization problem ([23]) can be interpreted as a refined 
condition of nonlinear dynamical stability (in that case 
5 is a Casimir or a pseudo entropy), as a sufficient condi- 
tion of thermodynamical stability in the MRS approach, 
or as a necessary and sufficient condition of thermody- 
namical stability in the EHT and NCD approaches (in 
these cases 5 is a generalized entropy). Noting that 5 
is proportional to the opposite of the potential enstro- 
phy T 2 = J q 2 dr, the maximization problem (|23|) is also 
equivalent to the phenomenological minimum potential 
enstrophy principle. In the following, to simplify the ter- 
minology, 5 will be called the "entropy" . 

The critical points of entropy at fixed circulation and 
energy are given by the variational principle 



SS - PSE - aST = 0, 



(28) 



where f3 and a are Lagrange multipliers that will be called 
"inverse temperature" and "chemical potential". This 
yields 



-/3ip — a. 



(29) 



We consider a domain of unit area A = 1 and we define 
(X) = J X dr. Then, we have a = — T — f3(ip) and we 
can rewrite the previous relation as 

q = -py,-{ i /,)) + r. (30) 

Substituting this relation in Eq. (j27j), we obtain 

-A^ + /3^ = r + /3(^)-/i, (31) 

with i\) = on the domain boundary. This is the funda- 
mental differential equation of the problem. It has the 
form of an inhomogeneous Helmholtz equation. 

Using Eq. ([30]) . the energy and the entropy can be 
written 

E = -\P («> 2 > - W 2 ) + \t{iP) - !(/#>, (32) 



\(? «^ 2 } - w 2 ) - ^r 2 . 



(33) 



Using Eq. ([32]) , an equivalent expression of the entropy 
is 



S = pE- l -pT^)+ l -p(W)- 1 - 



r 2 . 



(34) 



The last term in Eqs. ([33]) and ([34]) will be ignored in 
the following since it is just an unimportant additional 
constant (for fixed Y). 



B. The solution of the differential equation and the 
equation of state 

To study the maximization problem ([23]) . we shall fol- 
low the general methodology developed by Chavanis & 
Sommeria [22j. The main novelty with respect to their 
study is the presence of the topography h. We shall write 
the topography in the form h(x,y) = bH(x,y) where H 
is dimensionless. 

To study the differential equation 



A^ + pi> = r + pty) - bH, 



(35) 



we first assume that Y + P(ijj) ^ (i.e. ft/0) and we 
define 



i> 



r + W 



(36) 



We note that c = — b/a so that c plays the role of the 
inverse of the chemical potential (see Sec. IIIIC)) . There- 
fore, the condition a ^ is equivalent to c finite. With 
these notations, the differential equation ([35]) becomes 



A0 + fief) = 1 - cH, 



(38) 



with <j> = on the domain boundary. The solution of this 
equation can be written 



<p = <pl +C(p 2 , 

where <f>\ and <p2 are the solutions of 
-A0i+^i = 1, 

- A0 2 + fifa = -H, 



(39) 



(40) 



(41) 



with (f)i = and ^ = on the domain boundary. For 
given /?, the functions <pi and (\>2 can be obtained by solv- 
ing the differential equation numerically or by decompos- 
ing the solutions on the modes of the Laplacian operator 
(cf Appendix IC]) . For the moment, we assume that /3 is 
not equal to an eigenvalue j3 n of the Laplacian so that 
the solutions of Eqs. ([40]) and ([4T]) are unique and fi- 
nite. The cases /3 = j3 n must be studied specifically (see 
following sections). 

Taking the average of Eq. ([36]) and solving for (^), we 
obtain 



w = 



r(</>> 



Using Eqs. 



i - p(<t>) ■ 

and (03), the solution of Eq. (J3S]) is 



(42) 



ip 



Y4> 



i-W)' 



(43) 



Using Eqs. ([37]) and ([^2]) . the constant c is given by 



c= f (l-/3(</>», 



(44) 



where ((/>) itself depends on c. Substituting Eq. ([39 
Eq. ([44]) . we finally obtain 



in 



c = 



i-Wi> 



(45) 



and 



For a given normalized circulation Y/b, this relation com- 
pletely determines c as a function of j3. Therefore, the 
solution of Eq. ([35]) is given by Eq. ([43]) where is deter- 
mined by Eqs. ([39]) and ([45]) . We have thus completely 
solved the differential equation ([35]) for f3 ^ f3 n . We must 
now relate /3 to the energy. Substituting Eq. ([43]) in the 
energy constraint ([32]) . we get 



r + m ■ 



(37) 



(l-/3(0)) 2 = — ((0)-W 2 > 



c<0ff)). (46) 



Similarly, the entropy ([33]) can be written 

p 2 



25 

f2 



(1-W» : 



r((f)-(r)- 



(47) 



In the absence of topography (b = c = 0), we recover 
the equations of Chavanis & Sommeria [22J. In that 
case, there is a single control parameter A = T/V2E. 
In the present case, there are two control parameters: 
A = T/y/2E and \i = T/b. In order to have a well- 
defined limit T — )• 0, it is more convenient to take 
£ = 2E/b 2 = /i 2 /A 2 and /i = T/b as independent control 
parameters. Using Eq. (|55]h we find that the equations 
of the problem are given in a very compact form by 



2E 



1 



m-m 2 )-c{^H)\ 



%--?«* 



w 2 ), 



(48) 



(49) 



where (j) is given by Eq. ([39]) and c by Eq. ([45]) . 

Note that the r.h.s. of Eqs. ([^8]) and ([^9]) are functions 
of /? which can be easily computed numerically. There 
are two control parameters: the energy £ = 2E/b 2 and 
the circulation \i = T/b (normalized by the /3-effect pa- 
rameter b or by the amplitude of the topography). For 
the sake of simplicity we will denote, in the figures and 
in the discussion, E and T the energy and the circula- 
tion thus normalized (while b will be explicitly written in 
the formulae). For given T, Eq. ([48]) determines /3 as a 
function of E, i.e. the caloric curve /3(E). Of course, it 
is easier to proceed the other way round. We first fix T. 
Then, for each f3 we can determine c by Eq. ([^5]) and E 
by Eq. ([48]) to obtain E(f3). Inverting this relation we 
get /3(E) for fixed T. 

We now consider the case T + /3{i/j) = (i.e. a = 0). 
Equation ([35]) then becomes 



AiP + /3i/j = -bH, 



and the solution is 



v> = 



(50) 



(51) 



Substituting this relation in Eqs. ([32]) and ([33]) and using 
T + /3b (fa) = 0, we find that the energy and the entropy 
are given by 



2E 



-/3(cf> 2 2 ) - (faH) : 



(52) 



results of Chavanis & Sommeria [22j. This corresponds 
to a limit of large energies 1/E — >• 0. We expect geom- 
etry induced phase transitions between monopoles and 
dipoles. On the other hand, for f3 — » -foe, we recover 
the results of Fofonoff [4] . This corresponds to a limit of 
small energies 1/E — » -foe. In that case, the equilibrium 
state is unique and corresponds to a westward jet (in an 
antisymmetric domain with /^-effect H = y). 



C. The chemical potential 



The chemical potential is defined by 



a = -/3(V)-r. 



If /3(ip) + T = 0, then 



a = 0. 



If /3(tp) +1^0, according to Eq. (JS3), we have 



a 
~b 



(54) 



(55) 



(56) 



Therefore, up to a normalization constant, the chemical 
potential is equal to — 1/c. This gives to the parameter 
c a clear physical meaning. 

For a given value of the energy E, we can obtain the 
chemical potential curve a(T) in parametric form in the 
following manner. Fixing £", we find from Eqs. ([48]) and 
that c is related to j3 by a second degree equation 



f + /W2> 



-{{<hH)-{<h) + m<h<h))c + ^) 



k)= 0.(57) 



This determines c = c(E,/3). On the other hand, ac- 
cording to Eq. (pT5]h the circulation T/6 is related to /? 

by 



1 



c(S,/3) 



(1 - Pifa)) - P(<h 



(58) 



This determines T = T(E,/3). Therefore, for given E, 
these equations allow us to obtain c as a function of r/6 in 
parametric form with running parameter /?. This yields 
the chemical potential curve a(T). Then, we have to 
account for the particular cases where j3 is an eigenvalue 
f3 n of the Laplacian (see below). 



2S_ 
b 2 



-P 2 ((ti)-(fc) 2 ). 



(53) 



We can check that these equations are limit cases of the 
general equations ([48]) . ([49]) and ([39]) when c — » +oo. 
Indeed, the condition a —> is equivalent to c — >• +oo. 

Particular limits: It is interesting to mention the con- 
nection with previous works. For c — » 0, we recover the 



D. A critical circulation 

According to Eq. ([45]) . we note that the expression 
of c (related to the chemical potential a), involves the 
important function: 



^) = Wi)-i. 



(59) 



We shall call /3i k) the solutions of F(f3 {k) ) = and we 

shall denote simply /?*=/?* the largest solution. The 
function F(/3) and the inverse temperature f3* were in- 
troduced by Chavanis & Sommeria [22j . We will see that 
the temperature /3* plays an important role in the prob- 
lem. For f3 = /3*, we find that c = except if T = T* 
where T* is a critical circulation given by J42J: 



-y = -P*(<fo)*. 



(60) 



We will have to distinguish the cases T = T* and r^r*. 



E. The program 

We shall successively consider the case of an antisym- 
metric and a non-symmetric topography. It will be shown 
that the mathematical expressions simplify greatly for an 
antisymmetric topography so that it is natural to treat 
this case first. To be specific, we will consider a rect- 
angular domain and a topography of the form h = by. 
This corresponds to the situation studied by Fofonoff [4] 
in his seminal paper. It should therefore be given partic- 
ular attention. Then, we will consider a non-symmetric 
topography of the form h = b(y — yo) in a rectangular 
domain. Finally, it will be shown in Sees. IVl and IVTl how 
the results can be generalized to an arbitrary domain and 
an arbitrary topography. 

The inverse temperature ft is the Lagrange multiplier 
associated with the conservation of the energy E and the 
chemical potential a is the Lagrange multiplier associated 
with the conservation of the circulation Y. Thus, j3 = 
(dS/dE) r and a = (dS/dY) E . We shall first study the 
caloric curve f3(E) for a given value of the circulation T, 
then the chemical potential ot(Y) for a given value of the 
energy E. 



IV. THE CASE OF AN ANTISYMMETRIC 
LINEAR TOPOGRAPHY (FOFONOFF CASE) 

We consider a complete oceanic basin as in the study of 
Fofonoff Q. The domain is rectangular with — y^r/2 < 
x < y/r/2 and -l/(2y/r) < y < l/(2y/r) where r = 
L x /L y is the aspect ratio. The topography h = by (i.e. 
H = y) is linear and antisymmetric with respect to y = 0. 
This linear topography can also represent the /^-effect. 
The eigenvalues of the Laplacian in a rectangular domain 
will be called j3 mn (see Appendix IC]). Assuming (3 ^ /3 mri , 
it is easy to show from considerations of symmetry that 
02 is odd and 0i is even with respect to y. Therefore, 



2> = <0iy> = <0102> = 0. 



(61) 



When r + f3(ip) / (a / 0, c finite), the equations of 
the problem become 



ip= -(f>. 



2E _ 1 



c= f (1-/W), 



({<l>i)-P(<fi))-P{<l%)-(<hv), 



^ = -§Wi>-^i> 2 )-/3m 



(63) 



(64) 



(65) 



When r + f3(ip) = (a = 0, c — > oo), the equations of 
the problem become 



IE 



i> = b(p2 



-ml) - {<hv), 



25 
b 2 



^---P 2 (<t>i)- 



(66) 



(67) 



(68) 



01 + C0 2 ) 



(62) 



Finally, since (02 ) = 0, we find that r* = 0. We thus 
need to distinguish two cases depending on whether T = 
or T ^ 0. 



A. The caloric curve /3(E) for T = 

We shall first discuss the caloric curve /3(E) for T = 
r* = 0. Details on the construction of this curve can 
be found in Appendix IF II Since the multiple solutions 
occur for large values of E, it appears more convenient to 
plot ft as a function of 1/E like in the study of Chavanis 
& Sommeria [22]. The corresponding curve (3(E) can be 
deduced easily. 

For small energies, there exists only one solution of Eq. 
(|35|) and it is a global entropy maximum at fixed circu- 
lation and energy. On the other hand, for large energies, 
there exists an infinite number of solutions of Eq. (l35|) . 
i.e. there exists an infinite number of critical points of 
entropy at fixed circulation and energy. In order to select 
the most probable structure, we need to compare their 
entropies. For 1/E = we have S/E = j3 so we just 
need to compare their inverse temperature (3. More gen- 
erally, it can be shown that the entropy is a monotonic 
function of (3 (for given E and T) so, in the following, we 
shall select the solution with the largest f3. For large en- 
ergies, there is a competition between the solution with 
inverse temperature /?*, the solution with inverse tem- 
perature ^21 (where /?2i is the largest eigenvalue with 
zero average value (^21) = and (ip2iy) — 0) and the 
solution with inverse temperature /?i2 (where /3i2 is the 
largest eigenvalue with zero average value (^12) = and 
('0122/) 7^ 0)- I n t ne region where the solutions are in com- 
petition, the solution with the highest entropy is the one 
with j3 = maxj/3*, /32i, ^12}- The selection depends on 
the aspect ratio of the domain. In a rectangular domain 
elongated in the horizontal direction (r > 1), #21 > /?i2- 
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In a rectangular domain elongated in the vertical direc- 
tion (r < 1), /?i2 > /?2i- O n the other hand, as shown by 
Chavanis & Sommeria [22j, there exists a critical aspect 
ratio r c = 1.12 such that: max{/3*,/?2i,/3i2} = /?* for 
1/t c < r < r c , max{^*,/?2i,/?i2} = /?2i for r > r c and 
max{/3*,/32i,/3i2} = /?i2 for r < l/r c . To describe the 
phase transitions, we must therefore consider these three 
cases successively. 
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FIG. 1: Relationship between /3 and 1/E in a square domain 
[-0.5,0.5] x [-0.5,0.5] for T = (l/r c < r = 1 < r c ). The 
potential vorticity g is shown for several values of f3. For 
small energies we get Fofonoff flows. For large energies and 
1/t c < t < t c , the maximum entropy state is the monopole. 
There exists a second order phase transition between Fofonoff 
flows and monopoles (see zoom in Fig. [2}- Note that this 
second order phase transition exists only for r = 0. 

• l/r c < r < r c , as in Figs. Q]and[2j For small energies 
there exists only one solution of Eq. ([35]) and it is a global 
entropy maximum. For E — » 0, leading to f3 — » -f oo, we 
recover the classical Fofonoff solution (see Fig. [3]). Far 
from the boundaries, the Laplacian term in Eq. ([35]) can 
be neglected, which leads to /3-0 = T + /3(^) — % and 
u = — (b//3)x, representing a westward jet of velocity 
ii = b//3. The eastward recirculation can be obtained 
from a boundary layer approximation [4( . For intermedi- 
ate energies with /? > 0, the flow involves two symmetric 
gyres: the gyre at y > has positive PV and the gyre 
at i/ < has negative PV. For intermediate energies and 
f3 < 0, the situation is reversed: the gyre at y > has 
negative PV and the gyre at 2/ < has positive PV. All 
these solutions, forming the upper branch of the main 
curve, will be called Fofonoff flows J43j- They will be 
labeled (FP) and (FN) respectively. On the other hand, 
for large energies, there is a competition between several 
solutions of Eq. ([35]) . When l/r c < r < r Cl the max- 
imum entropy state is the solution with fi — /3*. The 
solution with j3 = /?* and 1/E = is a pure monopole 
as in the study of Chavanis & Sommeria [22j. It can 
rotate in either direction since the monopole (MP) with 
positive PV at the center and the monopole (MN) with 
negative PV at the center have the same entropy. For 
/? = /?* and 1/E > the monopole is mixed with a Fo- 



fonoff flow. It will be called mixed monopole /Fofonoff 
flow. For a fixed value of E, there exists two different 
solutions depending on the sign of c (see Fig. HJ. The 
caloric curve /3(E) displays a second order phase transi- 
tion marked by the discontinuity of q^(E) at E = E*(r). 
In conclusion, for 1/E > 1/E*(r) we have Fofonoff flows 
(with a(E) = 0), for < 1/E < 1/E*(t) we have mixed 
monopole/Fofonoff flows (with ±a(E) 7^ 0) and for 
1/E = we have a pure monopole (with a(E) — >• ±00). 




FIG. 2: Zoom of Fig. \T\in the region of second order phase 
transition. 




FIG. 3: Potential vorticity q in a square domain [—0.5, 0.5] x 



[-0.5, 0.5] for T = and f3 
so-called Fofonoff flow. 



10 . This corresponds to the 



• r > r c = 1.12 (horizontally elongated domains), as 
in Fig [5] For small and moderate energies, the situation 
is similar to that described previously (Fofonoff flows). 
On the other hand, for large energies, the situation is 
different. When r > r Cl the maximum entropy state is 
the solution with (3 = /?2i- The solution with f3 = #21 
and 1/E = is a pure horizontal dipole as in the study of 
Chavanis & Sommeria [22|. It can rotate in either direc- 
tion since the dipole (DP) with positive PV on the left 
and the dipole (DN) with negative PV on the left have 
the same entropy. For j3 = {3<i\ and 1/E > the dipole 
is mixed with a Fofonoff flow. It will be called mixed 
horizontal dipole /Fofonoff flow. The caloric curve 13(E) 
displays a second order phase transition marked by the 
discontinuity of -^(E) at E = £^21 (T = 0,r). In conclu- 
sion, for 1/E > l/E 2 i(T = 0,t) we have Fofonoff flows 
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izontal dipoles in the previous cases. In conclusion, for 
1/E > we have FofonofT flows (with a(E) = 0) and for 
1/E — » we have a pure vertical dipole (with a(E) = 0). 



FIG. 4: Potential vorticity q in a square domain [—0.5, 0.5] x 
[-0.5, 0.5] for T = 0, /3 = f3i 1] and (a) c = 0.5, (b) c = -0.5, 
(c) c = 1 and (d) c— — 1. Increasing values from blue to red. 



(with a(E) = 0), for < 1/E < l/E 2 i(Y = 0, r) we have 
mixed horizontal dipole/Fofonoff flows (with a(E) = 
and ±x(E)) and for 1/E = we have a pure horizontal 
dipole (with a(E) = and x ~~ ^ ±oo). 
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FIG. 5: Relationship between f3 and 1/E in a rectangular 
domain [-l/>/2, 1/V%] x [-l/(2>/2), l/(2>/2)] of aspect ratio 
r = 2 > r c for r = 0. For small energies we get Fofonoff flows. 
For large energies and r > r c , the maximum entropy state 
is the horizontal dipole. There exists a second order phase 
transition between Fofonoff flows and horizontal dipoles. 

• r < 1/t c = 0.893 (vertically elongated domains), as 
in Fig. [6] For 1/E > we recover Fofonoff flows as 
discussed previously. On the other hand, for 1/E = 0, 
when r < l/r c , the maximum entropy state is the pure 
vertical dipole with f3 = /3\2. In that case, there is no 
phase transition: the Fofonoff flows continuously form a 
vertical dipole for 1/E = 0. This can be explained by the 
fact that the vertical dipole does not break the symmetry 
of Fofonoff flows contrary to the monopoles and the hor- 
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FIG. 6: Relationship between f5 and 1/E in a rectangular 
domain [-l/(2>/2), l/(2>/2)] x [-l/v 7 ^ l/>/2] of aspect ratio 
r = 1/2 < 1/tc for r = 0. For 1/E > we get Fofonoff flows. 
For 1/E = and r < l/r c , the maximum entropy state is the 
vertical dipole. There is no phase transition. 

In Fig. [7J we plot the phase diagram in the (r, E) 
plane for Y = 0. Concerning the curve j3{E) at Y = 0, 
there is a second order phase transition between Fofonoff 
flows and horizontal dipoles for r > r c , a second order 
phase transition between Fofonoff flows and monopoles 
for l/r c < r < r c and no phase transitions for r < l/r c 
(the passage from Fofonoff flows to vertical dipoles is reg- 
ular). 



B. The caloric curve /3(E) for Y / 

We describe here the caloric curve /3(E) for Y ^ 0. 
Details on the construction of this curve can be found in 
Appendix IF 21 Three cases must be considered. 

• l/r c <r < r c , as in Figs. El M and [TOJ For 1/E > 0, 
the maximum entropy state is an asymmetric Fofonoff 
flow and for 1/E = 0, the maximum entropy state is 
the monopole /?*. Since there is no plateau, the caloric 
curve /3(E) does not display any phase transition: for 
r^0, /3(E) and §§j(E) are continuous. This is different 
from the case Y — 0. In conclusion, for 1/E > we 
have asymmetric Fofonoff flows (with a(E) ^ 0) and for 
1/E = we have a pure monopole (with a(E) = oo) 
rotating in either direction. Interestingly, when r — >■ 
(see Fig. [TO]) , the main curve is more and more "pinched" 
near the point (1/E = 0, f3 = /?*). This is consistent with 
the formation of a plateau (second order phase transition) 
when T = 0. 

• r > r c (horizontally elongated domains), as in Figs. 
[TT1 and [121 For small energies, the maximum entropy 
state is an asymmetric Fofonoff flow. For 1/E = 0, 
the maximum entropy state is the horizontal dipole /?2i- 
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FIG. 9: Zoom of Fig. 



FIG. 7: Phase diagram for an antisymmetric topography in 
the (t, E) plane when r = 0. Second order phase transi- 
tions arise, from monopoles to Fofonoff flows on the curve 
E — E*(r), and from horizontal dipoles to Fofonoff flows on 
the curve E — E2i(r). For 1/E = 0, the monopoles and the 
dipoles can rotate in either direction (the solution is degener- 
ate). The line Eo(r), corresponding to f3 = 0, is not plotted 
because it occurs for too high values of 1/E. This line sepa- 
rates the Fofonoff flows (FP) for f3 > to the Fofonoff flows 
(FN) for p < 0. 
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FIG. 8: Relationship between /3 and 1/E in a square domain 
[-0.5,0.5] x [-0.5,0.5] with T = 1 (l/r c < r = 1 < r c ). For 
1/E > 0, we get asymmetric Fofonoff flows. For 1/E = and 
l/r c < r < r c , the maximum entropy state is the monopole. 
Note that there is no plateau at f3 = f3* contrary to the case 
r = 0. Thus, for r/0, there is no phase transition. 



For intermediate energies, the maximum entropy state 
is a mixed horizontal dipole/asymmetric Fofonoff solu- 
tion. These solutions form a plateau j3 = /?2i- In 
that case, the caloric curve /3(E) displays a second order 
phase transition marked by the discontinuity of §§j(E) 
at E = E 21 (T,t). In conclusion, for 1/E > l/E 21 (T,r) 
we have asymmetric Fofonoff flows (with a(E) ^ 0), 
for < 1/E < l/£ , 2i(r,r) we have mixed horizontal 
dipoles/asymmetric Fofonoff flows (with a(E) ^ and 
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FIG. 10: Relationship between /3 and 1/E in a square domain 
[-0.5,0.5] x [-0.5,0.5] with T = 0.01. For r -> 0, the curve 
is more and more pinched near /3* explaining the formation 
of a plateau for F = 0. 



±x(i£)) and for 1/E = we have a pure horizontal dipole 
(with a(E) 7^ and \ ~^ +°°) rotating in either direc- 
tion. 

• r < l/r c (vertically elongated domains), as in Figs. 
[T3l and [TH For 1/E > 0, the maximum entropy state 
is an asymmetric Fofonoff flow and for 1/E = 0, the 
maximum entropy state is the vertical dipole /?i2. Since 
there is no plateau, there is no phase transition in that 
case. In conclusion, for 1/E > we have asymmetric 
Fofonoff flows (with a(E) ^ 0) and for 1/E = we have 
a vertical dipole (with a(E) ^ 0). 

In Fig. [15j we plot the phase diagram in the (r, E) 
plane for different values of T ^ 0. Concerning the curve 
/3(E) at fixed r^0, there is a second order phase tran- 
sition between Fofonoff flows and horizontal dipoles for 
r > r c and no phase transitions for r < r c (the passage 
from Fofonoff flows to monopoles and vertical dipoles is 
regular). 
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FIG. 11: Relationship between /3 and 1/E in a rectangular 
domain [-1/V% l/>/2] x [-1/(2 V^), l/(2>/2)] of aspect ratio 
r — 2 > r c with r = 1. For small energies, we get asymmetric 
FofonofT flows. For large energies and r > r c , the maximum 
entropy state is the horizontal dipole. There exists a second 
order phase transition between FofonofT flows and horizontal 
dipoles. 
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FIG. 12: Relationship between /3 and 1/E in a rectangular 
domain [-1/V% l/>/2] x [-l/(2\/2), 1/(2^2)] of aspect ratio 
r = 2 with T = 0.01. 



FIG. 13: Relationship between (3 and 1/E in a rectangular 
domain [-1/(2^/2), l/(2>/2)] x [-1/^2, l/\/2] of aspect ratio 
r — 1/2 < l/r c with r = 1. For 1/E > 0, we get asymmetric 
FofonofT flows and for 1/E — and r < l/r c , the maximum 
entropy state is the vertical dipole. There is no phase transi- 
tion. 
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FIG. 14: Relationship between (3 and 1/E in a rectangular 
domain [-1/(2^/2), l/(2>/2)] x [-l/>/2, l/>/2] of aspect ratio 
r = 1/2 with T = 0.01. 



C. The chemical potential curve a(T) for fixed E 

In the previous sections, we have studied the caloric 
curve /3(E) for a fixed circulation T. We shall now study 
the chemical potential curve a(T) for a fixed energy E. 
The general equations determining the chemical potential 
are given by Eqs. ([57]) and ([58]) . For an antisymmetric 
topography, using Eq. ([6T]) . the term in factor of c in Eq. 
([57]) vanishes and the foregoing equations reduce to 



1 



b 2 



■ P(<&) + (<hv) 



b 2 



b c(E,/3) 



i)-m) 



(i-Wi»- 



(69) 



(70) 



In that case, it is very easy to obtain the chemical po- 
tential curve a(r) for fixed E, parameterized by j3. This 
curve is antisymmetric with respect to V — >• — T. As in 
Sees. II V Al and IIVB] three cases must be considered. 

• l/r c <r < r c , as in Figs. [TBI and [171 For T = and 
1/E < 1/E*(t), we are on the plateau /3 = /?* of mixed 
monopoles/Fofonoff flows. There exists two solutions for 
each energy that have two symmetric values of the chem- 
ical potential ±a(E) ^ (see Fig. [4]). Their values are 
given by Eqs. ([56]) and ([57]) replacing f3 by /3*. At the 
end of the plateau, for 1/E = 1/E*(t), we get ±a = 0. 
For T = and 1/E > 1/E*(t), we are on the main 
branch and the chemical potential is equal to a(E) = 
for each energy. For r / there is a unique solution 
a(E) for each energy. In conclusion, for 1/E < 1/E*(r), 
the chemical potential curve displays a first order phase 
transition at T = marked by the discontinuity of a(T). 
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FIG. 15: Phase diagram for an antisymmetric topography 
in the (r, E) plane for different values of T / 0. Second 
order phase transitions arise, from horizontal dipoles to Fo- 
fonoff flows, on the curve E — E2i(Y,r). For 1/E = 0, the 
monopoles and the dipoles can rotate in either direction (the 
solution is degenerate). The line Eo(t), corresponding to 
f3 = 0, is not plotted because it occurs for too high values 
of 1/E. This line separates the Fofonoff flows (FP) for f3 > 
to the Fofonoff flows (FN) for f3 < 0. 



This corresponds to the transition from the monopole 
(MP) for T = 0+ to the monopole (MN) for T = 0". 
For 1/E > 1/E*(t), the curve ot(Y) is continuous and 
differentiate so there is no phase transition. 
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FIG. 16: Chemical potential as a function of the circulation 
in a square domain (l/r c < r = 1 < r c ). For 1/E < 1/.E*(t), 
there is a first order phase transition between positive and 
negative monopoles (here 1/E = 2 and 1/E*(t = 1) ~ 4.56). 

• r > r c , as in Figs. [18] and [19] For any Y and 
1/E < l/£ , 2i(r,r), we are on the plateau f3 = f3 2 \ 
of mixed horizontal dipole/Fofonoff flows. There exists 
two solutions for each energy (with ±x(2£)), but they 
have the same value of chemical potential a(E). For 
1/E > l/£ , 2i(r,r), we are on the main branch and 
the chemical potential takes a unique value a(E), corre- 
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FIG. 17: Chemical potential as a function of the circulation 
in a square domain (l/r c < r = 1 < r c ). For l/2£ > 1/E*(t), 
there is no phase transition (here 1/-E = 8 and 1/E*(t — 1) ~ 
4.56). 



sponding to a unique solution, for each energy. For T = 0, 
we always have a(E) = 0. Since £ , 2i(r,r) is a parabola 
of the form J5 2 i(r,r) = a(r)T 2 + c(r) with a > 0, its 
minimum value is E™ n (r) = £^21 (r = 0,r) obtained for 
T = 0. Let us now assume 1/E < 1/E^ n {r). Then, as 
long as |T| < T 2 i(E) (so that 1/E < l/E 2 i(T)), we are 
on the plateau j3 = /?2i and the chemical potential a(r) 
is a linear function of the circulation given by Eqs. ([56]) 
and ([58]) where /3 is replaced by /3 2 i- For |T| > r 2 i(£ ! ), we 
are on the main branch and ot(T) is given by Eqs. ([56]) . 
([57]) and ([55]). In conclusion, if 1/E < l/£;^ in (r), the 
chemical potential curve displays two second order phase 
transitions between horizontal dipoles and Fofonoff flows 
at T = ±T 2 i(E) marked by the discontinuity of ffr(r). 
If 1/E > l/£?2i m (T) there is no phase transition. 

• r < 1/t c , as in Fig. [20] For any T, we are on the 
main branch and the chemical potential takes a unique 
value a{E) for each energy. For Y = 0, we always have 
a(E) = 0. Furthermore the curve ot(Y) is continuous and 
different iable. In conclusion, there is no phase transition. 

It is interesting to mention the existence of particu- 
lar points in the phase diagram, (i) For l/r c < t < t Ci 
E = E* (t) is a critical point at which a first order phase 
transition appears (see Figs. [16] and [TTJ). (ii) For en- 
ergies such that 1/E < l/E 2 i(Y = 0,r c ) = 1/E*(t c ), 
(r = r c ,r = 0) is a bicritical point: for fixed energy, 
the system exhibits a bifurcation from a first order phase 
transition (see Fig. [T6]) to two second order phase tran- 
sitions (see Fig. [T8]) when increasing the aspect ratio of 
the domain, (iii) For r > r c , a second order azeotropy 
arises at E = £?2i(r = 0,r), where two second order 
phase transitions appear simultaneously from nothing 
(see Fig. [18] and [19]). The possible existence of these 
behaviors in systems with long range interactions was 
predicted by Bouchet & Barre [37j. It was evidenced 
by Venaille & Bouchet [23|, [24[ for Euler and geophysi- 
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FIG. 18: Chemical potential as a function of the circulation 
in a rectangular domain of aspect ratio r = 2 > r c . For 
1/E < l/EZl tn (r), there are two second order phase tran- 
sitions between horizontal dipoles and Fofonoff flows (here 
1/E = 500 and 1/E^ in (r = 2) « 1045). 
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FIG. 19: Chemical potential as a function of the circulation in 
a rectangular domain of aspect ratio r = 2 > r c . For 1/E > 
1/E^l in (r), there is no phase transition (here 1/E = 2000 and 
1/E^ in (r = 2) «1045). 



cal flows by extending the work of Chavanis & Sommeria 
j22j. We refer to these works for a more detailed de- 
scription of these phase transitions. We note that the 
critical point and the second order azeotropy are specific 
to flows with topography while the bicritical point also 
exists when h = 0. 



V. THE CASE OF A NON- SYMMETRIC 
LINEAR TOPOGRAPHY 

In the previous section, we have considered the case of 
a linear topography h = by that is antisymmetric with 
respect to the middle axis y = of a rectangular do- 
main (these results remain valid for more general anti- 



FIG. 20: Chemical potential as a function of the circulation in 
a rectangular domain of aspect ratio r — 1/2 < l/r c . There 
is no phase transition (here 1/E = 50). 



symmetric topographies of the form h = bH(y)). We now 
consider the case of a linear topography h = b(y — 2/o)? 
i.e. H = y — yo, with i/o / that is non-symmetric 
(the following results remain valid for more general non- 
symmetric topographies of the form h = bH(y)). We 
shall see that the details of calculations are a bit differ- 
ent while the structure of the main curves remains finally 
unchanged. 

For a non-symmetric topography, there is no particular 
simplification of the equations of the problem. Therefore, 
we must use the general equations of Sec. IIIIBI It has 
been shown in Sec. IIIIBI that there exists a critical circu- 
lation T*. This critical circulation depends on the aspect 
ratio r of the domain j44j and on the form of the topogra- 
phy. We shall consider successively the cases T ^ T*(t) 
and T = r*(r). 



A. The caloric curve /3(E) for T / T* 

The series of equilibria for a non-symmetric topogra- 
phy when r^T^ are similar to the ones obtained in Sec. 
IIVBI for an antisymmetric topography when r / (re- 
call that r* = in an antisymmetric domain). Details 
on their construction are given in Appendix IF 31 

• 1st case: l/r c < r < r c . For 1/E > we have asym- 
metric Fofonoff flows (with a(E) ^ 0) and for 1/E = 
we have a pure monopole (with a(E) = oo). The caloric 
curve /3(E) does not display any phase transition. 

• 2nd case: r > r c . For 1/E > l/E 2 i(T,r) we 
have asymmetric Fofonoff flows (with a(E) ^ 0), for 
< 1/E < I/E21 (r,r) we have mixed horizontal 
dipoles/Fofonoff flows (with a(E) ^ and Xi,2(E)) and 
for 1/E = we have a pure horizontal dipole (with 
a(E) 7^ and \ ~~ ^ +°°)- The caloric curve /3(E) displays 
a second order phase transition at the energy £?2i(I\t). 

• 3rd case: r < l/r c . For 1/E > 0, we have asym- 
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metric FofonofT flows (with a(E) ^ 0) and for 1/E = 
we have a vertical dipole (with a(E) ^ 0). There is no 
phase transition. 



B. The caloric curve /3(E) for r = T* 

The series of equilibria in a non- symmetric domain 
when r = T* are similar to the ones obtained in Sec. IIV Al 
for an antisymmetric domain when T = (recall that 
r* = in an antisymmetric domain). Details on their 
construction are given in Appendix IF 41 They can also be 
understood from the curves Y ^ T* of Sec. IV Al bv consid- 
ering the limit Y — >• T*. When Y — »• T*, the main curve 
is more and more "pinched" near the point (1/E = 0, 
f3 = /?*) and for Y = T* a plateau appears at tempera- 
ture /3 = /3* between 1/E = and 1/E = 1/E*(t). 

• If l/r c <r <r c : for 1/E > 1/E*(t) we have FofonofT 
flows (with a(E) ^ 0), for < 1/E < 1/E*(t) we have 
mixed monopoles/ FofonofT flows (with a\^(E) ^ 0) and 
for 1/E = we have a pure monopole (with a±^(E) —> 
±oo). The caloric curve /3(E) displays a second order 
phase transition at the energy E*(t). 

• If r > r c : for 1/E > l/E 2 i(Y*,r) we have FofonofT 
flows (with a(E) ^ 0), for < 1/E < l/E 2 i{T*,r) 
we have mixed horizontal dipoles/Fofonoff flows (with 
a(E) 7^ and xi,2(E)) and for 1/E = we have a pure 
horizontal dipole (with a(E) ^ and \ ~~ ^ +oo). The 
caloric curve /3(E) displays a second order phase transi- 
tion at the energy £ , 2i(r*,r). 

• If r < l/r c : for 1/E > we have FofonofT flows 
(with a(E) =/z 0) and for 1/E = we have a vertical 
dipole (with a(E) ^ 0). There is no phase transition. 



C. The chemical potential curve a(Y) for fixed E 

In a non-symmetric domain, the results are similar to 
those obtained in Sec. IIV CI even if the general equations 
are a little more complicated and the curve a(Y) is non- 
symmetric. 

• If l/r c < r < r c : for 1/E < 1/K,(r), there is a 
first order phase transition at Y = T* marked by the 
discontinuity of a(Y). For 1/E > 1/E*(t), the curve 
a(Y) is continuous and differentiate so there is no phase 
transition. 

• If r > r c : for 1/E < l/EV^ n (r), there are two 
second order phase transitions at Y — Y^AE) and 
T = r+(£) marked by the discontinuity of |g(r). If 
1/E > l/£^2i m (r), there is no phase transition. 

• If r < 1/t c : there is no phase transition. 

In conclusion: (i) for l/r c < r < r Cl E*(t) is a critical 
point; (ii) (r = r c ,r = Y*(t c )) is a bicritical point; (hi) a 
second order azeotropy arises for r > r c at E = £^ m (r). 



VI. SUMMARY AND GENERALIZATIONS 

We shall now generalize the previous results to the 
case of an arbitrary domain and an arbitrary topogra- 
phy h(x, y) = bH(x, y). Some illustrations of phase tran- 
sitions in geophysical flows with different topographies 
are given in [25j. We here develop the theory needed to 
interpret them. 

In the series of equilibria containing all the critical 
points of entropy at fixed circulation and energy, we must 
distinguish: 

• The main curve /3(E): each point of this curve cor- 
responds to a unique solution with a unique value of 
the chemical potential a. For 1/E — »• +oo (small en- 
ergies), leading to /3 — » -foo, we obtain FofonofT flows 
with an arbitrary topography. Far from the boundaries, 
we can neglect the Laplacian term in Eq. ([35]) leading 
to a stream function /3ip = Y + P(^} — bH and a velocity 
field u = 4z x ViJ. The recirculation at the boundary 
can be obtained from a boundary layer approximation. 

• The inverse temperature f3 = /?*: for Y ^ T*, there 
is only one solution with f3 = /3* that exists at 1/E = 0. 
This is a limit point of the main curve. For Y = T*, the 
solutions with (3 = (3* form at plateau going from 1/E = 
to 1/E*. On this plateau, each value of the energy 1/E 
determines two solutions with the same temperature /3* 
but with different chemical potentials ol\(E) and 01.2(E). 

• The inverse temperature f3 = f3'{ (largest eigenvalue 
with (ip") 7^ 0): this is a particular point of the main 
curve that corresponds to 1/E'{{Y) where E'/(Y) is a 
parabola. This point marks the domain of inequivalence 
between the grand canonical ensemble on the one hand 
and the canonical and microcanonical ensembles on the 
other hand (see (Hill and Sec. NTC\i . 

• The inverse temperature /3 = /3[ a (largest eigen- 
value with (ip[ a ) = and (Hi/j[ a ) = 0; eigenmode with 
zero mean orthogonal to the topography): these solutions 
form a plateau going from 1/E = to l/E[ a (Y) where 
E[ a (Y) is a parabola. On this plateau, each value of the 
energy 1/E determines two solutions with the same tem- 
perature /3[ a , the same chemical potential a(E) but with 
different mixing coefficients Xi(£") and X2(E). 

• The inverse temperature f3 = f3[ b (largest eigenvalue 
with (V4b) = and (Hi/j' lb ) ^ 0; eigenmode with zero 
mean non orthogonal to the topography): it exists only 
at 1/E = 0. This is a limit point of the main curve. 

The series of equilibria showing these different solu- 
tions is represented schematically in Fig. [21] for Y = Y* 
and in Fig. [22] for r / T*. In the domain of compe- 
tition, where there exists several solutions for the same 
energy, we must select the solution with the largest /3, 
which is the maximum entropy state. Thus, in this 
range, /3 = maxj/3*,/^,/^}. Note that the value of 
max{/3*, /3[ a , /3[ b } only depends on the geometry of the 
domain and on the form of the topography. 
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to /3n, /3[ a corresponds to /?2i and /3 16 corresponds to /3i2- 




FIG. 22: Relationship between £ and 1/E for r / T* 



A. The caloric curve /3(E) 



Let us first describe the caloric curve (3(E) for a given 
value of T. We need to distinguish three cases: 

1. If msix{f3*,/3[ a ,/3[ b } = (3*: (i) If T = r*, we have a 
second order phase transition since q§j(E) is discontin- 
uous at 1/E = 1/E*. (ii) If T ^ r*, there is no phase 
transition. 

2. If max{/3*,/3i a ,/?i 6 } = /^ a : we have a second order 
phase transition since §§}(E) is discontinuous at E[ a (T). 
This is because ^i a is orthogonal to the topography lead- 
ing to a symmetry breaking. However, this situation is 
not generic (see the remark at the end of lVIBp . 

3. If max{/3*, (3 r la , j3' lb ] = f3[ b : there is no phase transi- 
tion. This is because i(j[ b is not orthogonal to the topog- 
raphy, so that there is no symmetry breaking. 



B. The chemical potential curve a(r) 

Let us now describe the chemical potential curve a(T) 
for a given value of E. We need to distinguish three cases: 

1. If max{/3*, /3' la , (3[ b } = (3*: if 1/E < 1/E*, there 
is a first order phase transition at T = T* since a(T) is 
discontinuous at T = T*. For 1/E > 1/E*, there is no 
phase transition. 

2. If max{0„/3i o ,#J = p' la : if 1/E < l/{E[ a ) min , 
there are two second order phase transitions at V = 
T^(E) and T = r+(£). For 1/E > l/(E[ a ) min , there is 
no phase transition. 

3. If max{/3*, f3' la , /3[ b } = f3[ b : there is no phase tran- 
sition. 

All these results are fully consistent with those ob- 
tained by Venaille & Bouchet [23|, [24j using a different 
theoretical treatment. 

Remark: in general, the topography (e.g. in the 
oceans) is very complex and is generically not orthogonal 
to an eigenmode of the Laplacian, so that ifj[ a does not 
exist (note also that zero mean eigenmodes ijj' n exist only 
if the domain has specific symmetries which is generally 
not the case in the oceans). In addition, in generic situa- 
tions, r^r^. Therefore, in typical caloric curves /3(E), 
there are no plateaus at j3* and f3' la (see, e.g., |25j). Con- 
sequently, the interesting phase transitions (second order, 
bicritical points, azeotropy,...) described in [23|, [24{ and 
in this paper do not generically exist in geophysical flows 
(with complex topographies). One notorious exception of 
physical interest is the case of a rectangular basin with 
a linear topography (treated explicitly in Sees. [IV] and 
|V]) corresponding to the /^-effect, or the case of topogra- 
phies of the form h = bH(x) or h = bH(y). However, for 
more complex topographies there is no phase transition 
in the strict sense. Nevertheless, there always exists at 
least a "smooth" transition from a monopole to a dipole 
when we stretch the domain (at fixed high energy), and 
a "smooth" transition from a monopole or a dipole to a 
Fofonoff flow when we lower the energy (at fixed domain 
shape). On the other hand, (T = T*,r = r c ) is always at 
least a critical point at which a first order phase transi- 
tion appears (for sufficiently high energies). 



C. Thermodynamical stability 

In the previous sections, when several solutions were 
in competition for the same values of E and T, we have 
compared their entropies to select the maximum entropy 
state. It turns out that the maximum entropy state is 
the solution with the highest inverse temperature |22| . 
Therefore, if we consider fully stable states (global en- 
tropy maxima at fixed energy and circulation), the strict 
caloric curve corresponds to ft > max{/3*,/^ a , (3[ b }. We 
can also study the thermodynamical stability of the so- 
lutions by determining whether they are (local) entropy 
maxima or saddle points of entropy. Complementary sta- 
bility results have been obtained by Chavanis & Somme- 
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ria [22|, Venaille & Bouchet [23|,[2J and Naso et al. J2l|. 
We shall briefly recall their results and refer to the cor- 
responding papers for more details. 

Chavanis & Sommeria [22| and Naso et al. [2l| have 
obtained sufficient conditions of microcanonical instabil- 
ity by considering the effect of "dangerous" perturbations 
on the equilibrium states. Using their methods, it can be 
shown that the solutions with /3 < f3' la are unstable in 
the microcanonical ensemble. On the other hand, when 
fi'ia < A* ? ft can ^ e shown that the solution ijj' la corre- 
sponding to /? = f3' la and 1/E = is unstable in the 
microcanonical ensemble. By continuity, all the plateau 
j3 = f3[ a should be unstable (since the two extremities of 
this plateau are unstable). 

Venaille & Bouchet [23|, [24[ have shown that the solu- 
tions with j3 > max{/3*, j3' la1 f3' lh \ are stable in the canon- 
ical ensemble (hence in the microcanonical ensemble) 
while the other solutions are unstable in the canonical 
ensemble. However, as discussed by Chavanis & Som- 
meria [22[ and Naso et al. [2l| , some of these solutions 
can be metastable (local entropy maxima) in the micro- 
canonical ensemble. 

For /?>/?", the solutions are stable in the grand canon- 
ical ensemble (hence in the canonical and microcanonical 
ensembles) [23, 24]. This is related to the Arnold theo- 
rem (see, e.g., [lj). For max{/?*,/?i a ,/?iJ < /3 < #', 
the solutions are stable in the canonical (hence micro- 
canonical) ensemble but not in the grand microcanonical 
ensemble. This corresponds to a situation of ensemble 
inequivalence [23, [24J . 



As expected with such a low value of the energy (see 
Fig. [1]), the relaxation equation converges to a Fofonoff 
state. The inverse temperature j3 and the enstrophy are 
plotted as a function of time in Fig. [231 As expected, 
T2 monotonically decreases in time (equivalently, the en- 
tropy increases). On the other hand, the inverse tem- 
perature monotonically increases during the simulation. 
Both quantities remain constant once the steady state 
(Fofonoff flow, maximum of entropy) has been reached. 



*§■! 





FIG. 23: Potential vorticity at t = 0, 15, 25, 40, for T = 
and 1/E ~ 1050 (low energy) in a square domain. In that 
case, we obtain two rolls that are strongly influenced by the 
topography. For even lower energies, we obtain strict Fofonoff 
flows exhibiting a westward jet. 



VII. NUMERICAL SIMULATIONS 

We shall now perform numerical simulations to illus- 
trate the phase transitions described in the previous sec- 
tions. To that purpose, we use the relaxation equations of 
Sec. Ill CI that can serve as numerical algorithms to com- 
pute maximum entropy states with relevant constraints. 
We shall first integrate numerically Eq. (|20]h with the 
constraints (|21l22p . in an antisymmetric square domain 
with a linear topography and two different initial condi- 
tions such that T = 0. We use as boundary conditions: 



^\dD =0, 

q\dD = -a(i), 



(71) 
(72) 



where dD is the domain boundary. The first condition 
enforces free-slip on the boundary, while the second one 
is necessary for consistency of the steady state, charac- 
terized by Eq. ([29]) . 

We first integrate the relaxation equations with an ini- 
tial condition, q(x,y,t = 0), written as the sum of sine 
functions with random amplitudes and wave numbers 
ranging from 1 to 9 (see Fig. l23l a)). With such a field, 
T = and the energy is rather low: 1/E « 1050. The 
coefficient D is set to 0.3. The resulting density of po- 
tential vorticity at different times is plotted in Fig. [23j 



These results can be compared with those of Wang & 
Vallis [9]. In this study, the authors integrate the quasi- 
geostrophic equation with linear topography in an anti- 
symmetric square domain, with initial conditions (ran- 
dom eddies) and boundary conditions (free-slip) similar 
to ours. After about 10 eddy turnover times, the time 
averaged flow converges to a state close to the Fofonoff 
solution. Comparing our Fig. [23] with Fig. 4 of [9], it is 
clear that the dynamical behavior of the coarse-grained 
potential vorticity q solution of the relaxation equation 
is reminiscent of that of the time averaged q solution of 
the quasigeostrophic equation. In both cases, the eddies 
are first pushed to the western boundary, then two gyres 
form. These structures grow and fill out the northern 
and southern parts of the domain. Therefore, even if the 
relaxation equations are not supposed to provide a realis- 
tic parametrization of 2D turbulence, they may however 
give an idea of the true evolution of the flow towards sta- 
tistical equilibrium. It is worth noticing that, with the 
relaxation equations, the steady state has been reached 
after about 8000 time steps, while it took more than sev- 
eral ten thousands time steps for the quasigeostrophic 
equation to reach the statistically steady state. Indeed, 
by construction, the relaxation equations "push" the sys- 
tem in the direction of the statistical equilibrium state. 

We then impose an initial condition of high energy, 
1/E ~ 1.3.10 -5 (see Fig. I25la)) and integrate numeri- 
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FIG. 24: Time evolutions of the inverse temperature (3 and 
of the enstrophy T2, for r = and 1/E c± 1050 (low energy) 
in a square domain, corresponding to Fig. [23l 



cally Eq. (|20]h with the constraints ([21122ft . Since we are 
only interested in the final state, and not in the way the 
system converges towards it, we do not implement the ad- 
vection term. The coefficient D is set to 1. The q density 
is plotted at different times in Fig. [25j and the time evo- 
lutions of the inverse temperature and of the enstrophy 
are shown in Fig.[26j We find that the system first relaxes 
spontaneously towards a horizontal dipole (see Fig. I25f b) 
and first plateaus of Fig. [26]) . In the absence of external 
perturbation, the system remains in this state for a long 
time, even if it is predicted to be unstable (see Fig. [2j). It 
can be noticed that the value of the inverse temperature 
is slightly smaller than /3i2, confirming that the system 
is on the lower branch of Fig. [2j As already shown in the 
case of the 2D Euler equations [21], we illustrate here that 
saddle points can be very long-lived. Inspired by the sta- 
bility analysis performed in [21[, we add to the system, 
at t = 10, a perturbation of the form Sq = 1 — /3*0i,*, 
where 0i } * is the solution of Eq. ([40]) for /3 = /?* (see 
|21|). As expected (see Fig. [2]), the system is then imme- 
diately destabilized, and relaxes towards the monopole 
which is the maximum entropy state (Fig. l25l c-d) and 
second plateau of Fig. [26]). Depending on the sign of the 
perturbation, the system can relax towards the direct or 
towards the inverse monopole. 

To summarize, we have illustrated the fact that in an 
antisymmetric square domain with linear topography, if 
T = 0, the stable state is strongly correlated to the to- 
pography at low energy (Figs. [23] and [24]) . but is not 
influenced by it at high energy (Figs. [25] and [26]) . At 
high energy, it is influenced by the domain geometry: in a 
square domain, we get a monopole whereas in a rectangle 
sufficiently elongated in the x direction (r > r c = 1.12), 
we get a dipole. To illustrate this remark, we integrate 
numerically Eq. ([20]) , with the constraints ([21l22p , in an 
antisymmetric rectangular domain of aspect ratio r = 2, 
starting from an initial condition of zero circulation and 
of high energy (1/E ~ 1.6.10 -5 ). The resulting potential 
vorticity density is plotted at different times in Fig. [27j 
The system first converges towards the vertical dipole 
(saddle point), with f3 ~ f3\2 (see Fig. l27T b) and first 
plateaus of Fig. [28]) . It then destabilizes spontaneously, 
and converges towards the horizontal dipole (stable state) 
(see Fig. l2?Td) and second plateaus of Fig. [28]), with 
(3 ~ /?2i- These results can be compared to Fig. [5] 
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FIG. 25: Potential vorticity at t = 0, 10, 90, 140, for F = 
and 1/E ~ 1.3.10 -5 (high energy) in a square domain. The 
system first relaxes towards the dipole (saddle point). At 
t = 10, the optimal perturbation is applied to the system. It 
then converges towards the monopole (stable steady state). In 
that case, the stable state is not influenced by the topography. 




FIG. 26: Time evolutions of the inverse temperature /3 and 
of the enstrophy r 2 , for T = and 1/E ~ 1.3.10" 5 (high 
energy) in a square domain, corresponding to Fig. [25j At 
t = 10 (indicated by the arrows), the optimal perturbation is 
applied to the system. 



Other numerical simulations of the relaxation equa- 
tions illustrating phase transitions in geophysical flows 
with nonlinear topography are reported in [25j . Interest- 
ingly, the initial conditions of high energy that we have 
considered in the present paper are similar to those used 
in Fig. 7(c) (square domain) and 8 (rectangular domain 
of aspect ratio r = 2) of [25]: q(t = 0) is proportional 
to ^44 and the energies are identical. We find that the 
final states are the same with both topographies. This 
is to be expected since the topography should not influ- 
ence the maximum entropy state in the limit of high en- 
ergy. However, we observe that the relaxation equations 
with the nonlinear topography of [25] directly converge 
towards the equilibrium state whereas, with the linear 
topography, the system first converges towards the verti- 
cal dipole (saddle point). Therefore, at very high energy, 
while the topography does not influence the stable states, 
it seems to play a role in the dynamics of the relaxation. 
Of course, this conclusion is reached on the basis of our 
relaxation equations. It would be interesting to know 
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0, 10, 21, 60, for T = 



5 (high energy) in a rectangular domain of 
The system first relaxes towards the ver- 



FIG. 27: Potential vorticity at t 

and 1/E~ 1.6.10"" 

aspect ratio r = 2 

tical dipole (saddle point). It then destabilizes spontaneously 

and converges towards the horizontal dipole (stable steady 

state). In that case, the stable state is not influenced by the 

topography. 
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raphy and only weakly by the domain geometry. For 
large energies, we obtain geometry induced phase tran- 
sitions between monopoles and dipoles like in the study 
of Chavanis & Sommeria [22J. In that case, the flow 
is strongly influenced by the domain geometry and only 
weakly by the topography. In rectangular domains elon- 
gated in the x direction (r > r c ), as we decrease the 
energy of the flow, we describe symmetry breaking phase 
transitions between horizontal dipoles and Fofonoff flows. 
Alternatively, in rectangular domains elongated in the y 
direction (r < l/r c ), the smooth transitions between ver- 
tical dipoles and Fofonoff flows occur without symmetry 
breaking. This phenomenology, illustrated in a rectangu- 
lar domain with a linear topography, has been generalized 
to arbitrary domains and arbitrary topography. 

Our study returns and confirms the results previously 
obtained by Venaille & Bouchet [23|, [24[ by a different 
method. These authors provide a very detailed statisti- 
cal analysis of the problem, emphasizing the notions of 
bicritical points, azeotropy and ensemble inequivalence. 
Our approach, that is less abstract and illustrated by 
several explicit calculations, provides a useful comple- 
ment to their study. It gives another way of describing 
the complicated and rich bifurcations that occur in geo- 
physical flows. Our theoretical results extend the work 
of Chavanis & Sommeria [22| to the case of geophysi- 
cal flows (with a topography) and allow to interpret the 
phase transitions studied numerically by Chavanis et al. 
[251 ] in the case of complex topographies. 



FIG. 28: Time evolutions of the inverse temperature /3 and of 
the enstrophy r 2 , for T = and 1/E ~ 1.6.10" 5 (high energy) 
in a rectangular domain of aspect ratio r = 2, corresponding 
to Fig. [271 



whether it remains valid for more realistic parameteriza- 
tions. 



VIII. CONCLUSION 

In this paper, we have studied the nature of phase tran- 
sitions in simple models of oceanic circulation described 
by the QG equations with an arbitrary topography h. 
We have assumed a linear relationship between potential 
vorticity and stream function corresponding to minimum 
potential enstrophy states. We have given several inter- 
pretations of this minimum potential enstrophy principle 
in connection with statistical mechanics, phenomenolog- 
ical selective decay principles, and nonlinear dynamical 
stability. We have explicitly treated the case of a rect- 
angular basin and an antisymmetric linear topography 
h = by like in Fofonoff 's classical paper. For small en- 
ergies, we recover Fofonoff 's westward jet solution. In 
that case, the flow is strongly influenced by the topog- 



Appendix A: Minimum potential enstrophy principle 

We shall here discuss the difficulty to justify a mini- 
mum potential enstrophy principle based on the viscosity. 
In the presence of viscosity, the QG equations become 



dq 

m 



u • Vg = z/Acj, 



(Al) 



where q = uj + h. We emphasize that the quantity dis- 
sipated by viscosity is the vorticity uj, not the potential 
vorticity q. Therefore, the rate of change of potential 
enstrophy Y2 = J q 2 dr is 



' / qAou dr = 2u 



T 2 = 2v I qAuj dr = 2v qA(q - h) dr. (A2) 

For the 2D Navier-Stokes equation (h = 0), we obtain 
after an integration by parts 

f 2 = -2i/ f(\/uj) 2 dr<0, (A3) 

so that the enstrophy decreases monotonically under the 
effect of viscosity. However, for the viscous QG equa- 
tions, we do not have a monotonic decay of potential 
enstrophy. Following Bretherton & Haidvogel [13|, we 
must assume h <C q in Eq. (|A2j) in order to have 
a monotonic decay of potential enstrophy. In more 
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general cases, the phenomenological minimum potential 
enstrophy principle is not clearly justified. These re- 
marks also apply to the functionals ([7]). Indeed, while 
these functionals increase monotonically under the ef- 
fect of viscosity for the 2D Navier- Stokes equations since 
S = v J C"(uj)(Vuj) 2 dr > 0, their monotonic increase is 
not guar anted for the QG equations unless ft <C q. 

Therefore it is difficult to justify a principle of the form 
(j6j) for the viscous QG equations. By contrast, the in- 
terpretations of this principle that we have given in the 
framework of the inviscid QG equations are valid. 



Einstein- type formula £ = Df3(t) involving the diffusion 
coefficient and the inverse temperature. At equilibrium, 



Au 



P 



C"(") 



u. 



(B6) 



This equation can be directly derived from the relation 
C (uo) = — /3ip — a determining the steady states of Eq. 
(jBip . Combining Au = z x Vuo, Vuo = oj' (ijj)V ijj and 
C"(uo) = — /3 J 'uo' '(i/j) , we recover Eq. (|B6|) . In particular, 
if S is the opposite of the enstrophy S = —(1/2) J uo 2 <ir, 
the relaxation equation (|B5|) becomes 



Appendix B: Link with the Neptune effect 

In this Appendix, we show that the relaxation equa- 
tions of Sec. Ill CI are consistent with the Neptune effect 
of Holloway [36|. This discussion is related to the one 
given by Kazantsev et al. [10] but the justification of the 
relaxation equations that we consider is different. 

Let us first consider the case ft = and R = oo and 
derive the relaxation equation for the velocity field. We 
start from the equation for the vorticity field (|9]) that 
becomes 



duo 
~dt 



u • Vu; = V • 



D\Vuo 



C"(oo) 



V^ 



(Bl) 



where we recall that D can depend on position and time. 
For a 2D field, we have the identity V x (z x a) = ( V • a)z. 
Therefore, we can rewrite the foregoing equation as 



dut 

di 



— + u • Vw 



Vx 



z x D Vw 






Vip 



(B2) 



The corresponding equation for the velocity field is 



| + (u.V)u = -Iv, + zxZ)(v W , 



VV> 



(B3) 



where p is the pressure and p the density. To pass from 
Eq. (pEJ) to Eq. ([B3]), we have used the identity (u • 

V)u = V(u 2 /2) -uxwz and the identity V x (u x ujz) = 
— (u • V)cjz that is valid for a 2D incompressible flow, so 
that finally V x [(u • V)u] = (u • V)gjz. Now, using 
u = — z x V^ and the identity 

Au = V(V • u) - V x (V x u) 

= -V x (ujz) = z x Vo;, (B4) 

valid for a 2D incompressible flow, we finally obtain 



^ + (u.V)u = -iv P+ ^(Au- 



(B5) 



We see that the drift term in the equation for the vorticity 
takes the form of a friction in the equation for the veloc- 
ity. Furthermore, the friction coefficient is given by an 



<9 U , ^ 
w + (u.V)u 



— Vp + I>(Au-0(t)u). (B7) 
P 



At equilibrium, Au = /3u. 

Let us now consider the relaxation equation (|2j) includ- 
ing a topography h. It can be rewritten 



doj 
~dt 



u • Vq = V • 



D Vw + Vft • 



C"(q) 



V^ 



• (B8) 



Proceeding as before, the corresponding equation for the 
velocity field is 



- + (u.V)u- 



hz x u 



— Vp- 
9 



D Au + zxV/i 



C"(q) 



(B9) 



where we have used V x (z x h\i) = (V • (hu))z = 
(u • Vft)z for a 2D incompressible velocity field. In par- 
ticular, if S is the opposite of the potential enstrophy 
S = —(1/2) J q 2 dr, the relaxation equations (|B8|) and 
(lB9l) become 



u • Vq = V • [D (Vcj + Vft + f3(t)V<i/j)] , (BIO) 



duu 



and 

-- + (u-V)u + ftzxu = — Vp4 
at p 

D (Au + z x Vft-/3(£)u). 

Inspired by the original idea of Holloway [36| , 
duce a velocity field based on the topography 



(Bll) 

we intro- 



1 



m 



z x Vh, ip* 



h 

wr 



m^ 



With these notations, the equation for the velocity field 
can be rewritten 



du , ^. 7 1 

— + (u • V)u + ftz x u = —Vp- 
ot p 

D Au-D/3(£)(u-u*). 



(B13) 



It involves a turbulent viscosity and a friction force pro- 
portional to the difference between the velocity field 
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u and the velocity field u* based on the topography. 
The friction coefficient is given by an Einstein relation 
£ = Dfi. This term "pushes" the flow towards the topo- 
graphic flow u* . This corresponds to the so-called "Nep- 
tune effect" of Holloway [36]. On the other hand, the 
diffusion term allows some deviation with respect to the 
topographic flow. At equilibrium 



Au 



(u - u*). 



(B14) 



If D is constant, the equation (jB10|) for the vorticity can 
be written 

^ + u • Vq = DAuj - Dfi(t)(uj - cj # ). (B15) 



At equilibrium, 



Auj = fi(uj — oj*). 



(Bi6) 



If we neglect the Laplacian, we get uo = uo* . This is valid 
in the limit fi — » +00. This is equivalent to neglecting 
the Laplacian in the fundamental differential equation 
— At/j + h = —fitjj — a. This leads to if) = — i(/i + a), 
equivalent to u = u* and uo = uo*. As we have seen, 
this corresponds to the standard Fofonoff [4] flows that 
are completely determined by the topography (far from 
the boundaries). More generally, Eq. (|B16|) takes into 
account finite temperature effects that can induce devia- 
tions to the standard Fofonoff flows. These finite temper- 
ature effects (corresponding to sufficiently high energies) 
are precisely those that have been described in this paper. 
In conclusion, the relaxation equation (J9]) derived from 
the Maximum Entropy Production Principle (MEPP) 
|26j is relatively consistent with the oceanographic 
parametrization of Holloway [36], especially when the 
(generalized) entropy is the neg-enstrophy. This is in- 
teresting because the parametrization of Holloway J36J 
has been used in realistic oceanic modeling where the 
Neptune effect was shown to play a significant role. This 
suggests that our parametrization can be of relevance also 
in the physics of the oceans. 



Appendix C: Modal decomposition 

We define the eigenfunctions and eigenvalues of the 
Laplacian by 



At/Vi = Mn 



(CI) 



with t/Vi = on the domain boundary. These eigenfunc- 
tions are orthogonal and normalized such that (f/ViVVn) — 
S nm . Since — /(V^n) 2 ^ 1 * = fin J ^n^, we n °te that 
fin < 0. Following Chavanis & Sommeria [22j, we dis- 
tinguish two types of eigenmodes: the odd eigenmodes 
ip n such that (ip n ) = and the even eigenmodes ip!^ such 
that (ip'^ 7^ 0. We note fi n and fi" the corresponding 
eigenvalues. 



In a rectangular domain of unit area whose sides are 
denoted a = y/r and b = 1/y/r (where r = a/b is the 
aspect ratio), the eigenmodes and eigenvalues are 

Vw = 2 sin(m7r(x/vr + 1/2)) s'm(n7r(y/ry + 1/2)), (C2) 



2 f™ 2 2 

fimn = -7T I V TU 



(C3) 



where the origin of the Cartesian frame is taken at the 
center of the domain. The integer m > 1 gives the num- 
ber of vortices along the x-axis and n > 1 the number of 
vortices along the y-axis. We have (^ m n) — if m or n 
is even and (^mn) 7^ if m and n are odd. The largest 
eigenvalue with non zero mean is fin. The largest eigen- 
value with zero mean is ^21 for r > 1 and /?i2 for r < 1. 
We want to solve the differential equations 



and 



A0i+/30i = 1, 



A0 2 + fifc = -H, 



(C4) 



(C5) 



with (j)\ = and 02 = on the domain boundary. This 
can be done by decomposing the solutions on the eigen- 
modes of the Laplacian using / = ^ mn (fiimn}^mn- 

For fi ^ fimn, the solution of Eq. (|C4|) is unique and 
given by 



E 



(A 



From this relation, we obtain 



-Vvi 



1) 



E 



EyPmn/ 

P-Pr, 

mn 



0>r, 



(P " Pranf 



d(4>i) 
dp ' 



^> = E 



{^mn){H^mn) 
P Pmn 



(C6) 



(C7) 



(C8) 



(C9) 



For fi 7^ fimn, the solution of Eq. (jC5|) is unique and 
given by 



E 



(Hlpmn) 



Ymn- 



,J~Pn 

From this relation, we obtain 



!*> = "£ 



P-Pn 



mn 



(P - P mn ) 2 ' dp ' 



(CIO) 



(Cll) 



(C12) 
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(#Vw)(Vw) 
P Pmn 



We also have 



n<P2 



£ 



(Htp mn )(^ mn ) 
(P-Pmn) 2 ' 



Finally, we remark that 
(4>iH) 



(C13) 



(C14) 



(C15) 



On the other hand, if the topography is antisymmetric 
with respect to y = (as in the case of a linear topogra- 
phy H = y), we have (#Vw)(Vw) = implying 



hH) = (0 2 ) = (0102) = 0. 



(C16) 



In that case, 0i is even and 02 is odd with respect to the 
variable y. This can be directly seen at the level of the 
differential equations (jC4|) and (jC5|) . 



Appendix D: The case (3 = 

It can be interesting to consider the case (3 = specif- 
ically. This regime of infinite temperatures corresponds 
to uniform potential vorticity q = — a. The differential 
equation ([35|) reduces to 



- Ai/j = T - bH. 

The solution is -0 = T0i + 602 where 
solutions of Eqs. (gD]) and (gT} with /? : 
for the energy (|48|) can be written 



(Dl) 

!>i and 02 are the 
= 0. The equation 



2E 



iH)> 



b 2 H), (D2) 



where we have used c = b/Y according to Eq. ([^5]) . For 
given T, Eq. (|D2|) determines the energy E(Y) for which 
(3 = (i.e. g is uniform). It is given by a parabola. If 
the topography is antisymmetric with respect to y = 0, 
using Eq. (|C16|) , the foregoing equation reduces to 



2E u \ r ' 



iH). 



(D3) 



Appendix E: The low energy limit 

In the limit (3 — >• -foo, a boundary layer approxima- 
tion can be made. As a first approximation, the Lapla- 
cian can be neglected in the differential equations (40) 
and (41), everywhere in the domain except close to the 
boundary. To leading order, the function 0i is given by 
0i = 1/(3. The correction, close to the boundary, behaves 
like Ae~vP^, where £ is a coordinate perpendicular to 
the boundary pointing towards the inside of the domain. 
The constant A is determined by the condition 0i = on 



the boundary. Similarly, to leading order, the function 
02 is given by 02 = —y/(3 with a correction scaling like 
A e -VPC dose to the boundary. 
In this limit, we obtain 



<#> = 

Similarly, 

(022/) = 
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Finally, substituting these results in the energy equation 
(64), we obtain for (3 — >• +oo: 



2E 



1 



1 



1 



(E5) 



Appendix F: Technical construction of the caloric 
curves 

1. Antisymmetric linear topography and r = 

Several cases must be considered in order to construct 
the caloric curve. In the first one, V = and (ip) = 
(i.e. a = 0, c = oo). Since ip = 602 with (02) = 0, 
these solutions exist for any value of (3. The curve (3(E) 
relating their temperature to the energy is given by Eq. 
(f67|) . This forms the main curve (see Figs. [TJ [2j [5] and 
[6]) . The entropy of these solutions is given by Eq. (|68|) . 

Another possible situation is Y = and (ip) ^ (i.e. 
a 7^ 0, c finite). According to Eq. ([63]) . we see that c 

is finite iff 1 - /3(0i) = 0, i.e. (3 = fii k) . Therefore, the 
temperatures of these solutions take only discrete values. 
The value of c is determined by the energy E according 
to Eq. ([51]) with (3 replaced by (3* . This determines 
two solutions ±c(E) (i.e ±a(E)) for each value of the 
energy (see for instance Fig. 3]). The ensemble of these 

solutions form a plateau (3 = /3* (see Figs. [TJ [2j [5] and 
[6]) starting at 1/E = (c = 0, a = ±oo) and connecting 
the main curve at 1/E = l/E'i (r) (c = ±oo, a = 0). 
When 1/E = 0, we find S/E = (3i k) . On the other hand, 

£?* (t) is given by Eq. ([67j) with /3 replaced by /3* . 
The evolution of E*(r) with r is shown in Fig. 



23 



Finally, it will be assumed that f3 is equal to an eigen- 
value j3 mn of the Laplacian. In that case, (ip) is neces- 
sarily equal to (i.e. a = 0, c — » oo), which corresponds 
to ip — bfa. Four cases must be distinguished: 

• Case 1: /3 = j3 mn with (yipmn) ^ and (Vw) ^ 0. 
This is not possible for an antisymmetric topography. 

• Case 2: fi = j3 mn with (yi/jmn) ^ and (Vw) = 
(m odd and n even) like /?i2. In that case, Eq. ([4T]) has 
no solution for f3 = f3 mn and fa diverges like 



(yipmn) 
P Pmn 



Ipn 



(Fl) 



for f3 — »> /3 mn . Then, we find from Eqs. ([67]) and ([68]) that 
1/E" — » and S/E —> (3 mn . This is a limit case of the 
main curve. Note that the sign of fa (hence q) changes 
as we pass from P mn to /3~ n . 

• Case 3: p = /3 mn with (# mn ) = and (Vw) = 
(m even and n arbitrary) like /?2i- In that case, the 
solution of Eq. (|H|) with /? = /3 mn is not unique. To the 
solution (jCllj) of Eq. ([4T]) with /3 = /3 mn we can always 
superpose the eigenmode ^ mn with an amplitude x- This 
corresponds to a "mixed solution" 



fa -» 02 + X^n 



(F2) 



which has zero average as required. The amplitude x 1S 
determined by the energy constraint ([67]) with /3 replaced 
by Pmn- This determines two solutions ±x(E) for each 
value of the energy (note that these two distinct solutions 
have the same value of the chemical potential a = 0). 
The ensemble of these mixed solutions form a plateau 
P = Pmn (see Figs. [TJ H [5] and |6j) starting from 1/E = 
(x = =Loo) and connecting the main curve at 1/E = 
l/E mn (T = 0,r) (x = 0). For 1/E = 0, we find S/E = 
Pmn- On the other hand, E mn (T = 0,r) is given by Eq. 
([67]) with /? replaced by /? mn . The evolution of £"21 (r = 
0, r) with r is shown in Fig. [71 Note that in a square 
domain (r = 1), since Pmn = P nrn by symmetry, we find 
that 1/E = for m even and n odd according to Case 2 
(see Figs. [T]and[2j). In that case, the plateau reduces to 
a point. 

• Case 4: p = p mn with (yi/j m n) = and (Vw) ^ 
(m odd and n odd) like P\\. In that case, the solution 
of Eq. ([4T]) with /3 = /3 mn is not unique and we expect a 
mixed solution of the form fa ^ fa -\- X^mn- However, 
ip = bfa has a zero average value (as required) only for 
X = 0. Hence, we are just left with a limit case of the 
main curve. It corresponds to the energy E mn (T = 0,r) 
given by Eq. ([67]) with P replaced by P mn . 



2. Antisymmetric linear topography and T / 

Once more, several cases must be considered. It can 
be first assumed that r^0 and T-\-p(if;) = (i.e. a = 0, 
c — >> 00), which corresponds to ip = bfa. Since (fa) = 0, 
there are no such solutions for T 7^ and P 7^ Pmn- 

Another possibility is V 7^ and V + /3(ip) 7^ (i.e. 
a 7^ 0, c finite). These solutions exist for any value of 



p. The curve (3(E) relating their temperature to their 
energy is given by Eqs. ([63]) and ([6l]) . This forms the 
main curve (see Figs. I5BT4]). Their entropy is given by Eq. 
(165]) . For p = pi k \ we have c = leading to 1/E = 0. 
In that case, S/E — /3* . The caloric curve P(E) is 
unchanged when T — >• — T since only T 2 appears in the 
equation for the energy. However, the sign of c (hence of 
a) changes when r — » — T so that the flow structure is 
different. In the figures, we consider T > 0. 

Finally, if P is equal to an eigenvalue P mn of the Lapla- 
cian, four cases must be distinguished. 

• Case 1: p = p mn with (# mn ) 7^ and (Vw) ^ 0. 
This is not possible for an antisymmetric topography. 

• Case 2: P = p mn with (yipmn) ^ and (Vw) = 
(m odd and n even) like P12. We are necessarily in the 
case T + /3(ip) 7^ (i.e. a 7^ 0, c finite). We see that Eq. 
(jJTJ) has no solution for (3 = Pmn and that fa diverges 
like 



, (y^mn) , 

P Pmn 



(F3) 



when /3 — >> /3 mn . In that case 1/E — >> and 5/^ — > /3 mn - 
This is just a limit case of the main curve. Note that the 
sign of fa (hence q) changes as we pass from P mn to P mn . 
• Case 3: P = P mn with (yipmn) = and (Vw) = 
(m even and n arbitrary) like /?2i- We are necessarily in 
the case T + /3(ip) 7^ (i.e. a 7^ 0, c finite). The solutions 
of Eqs. ([40]) and (|4T]) are not unique since we can always 
superpose to Eqs. ([C6p and (|C10p an eigenfunction ^ mn 
of the Laplacian. Thus 



fa -> fa + X^n 



</>2 ->" 02 + X^n 



(F4) 
(F5) 



The amplitude \ i s determined by the energy constraint 
([M]) with /3 mn replacing /3. This determines two solutions 
±x(£') for each value of E (note that these two distinct 
solutions have the same value of c(E), hence of a(E), 
given by Eq. ([63]) ). The ensemble of these mixed solu- 
tions form a plateau /3 = Pmn (see Figs. [8lfH]) starting 
from 1/S = (x — >• +00) and reaching the main curve 
at 1/E = l/E mn (T,r) (x = 0). For 1/E = 0, we find 
S/E — » /3 mn . On the other hand, £^ mn (r,r) is given by 
Eqs. ([63]) and ([M]) replacing /3 by /3 mn . It corresponds to 
a parabola of the form E mn (T : r) = a(r)T 2 + c(r) with 
a > (since (0i) - /3(0?) > Q). The evolution of 
^2i(r,r) with r is shown in Fig. [15] for different val- 
ues of T. Note that in a square domain (r = 1), since 
Pmn = Pnm by symmetry, we find that 1/E = for m 
even and n odd according to Case 2 (see Figs. I51ITQ]). In 
that case, the plateau is reduced to a point. 

• Case 4: p = p mn with (yipmn) = and (Vw) ^ 
(m odd and n odd) like p\\. If we consider T + /3(^) 7^ 
(i.e. ft/0,c finite), we see that Eq. ([40]) has no solution 
for P = /3 mn and that fa diverges like 



, (Vw) , 

P Hmn 



(F6) 
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for /3 -)► /3 mn . 

-(Vr)/3(Vw>7(/5 

a finite value 



Then we find that c ~ 
Pmn) —> oo and that E tends to 



2E 



1 



mmn) 2 & 



/?<02> - (022/>, 



(F7) 



with /3 = /3 mn . It corresponds to a parabola of the form 
^ mn (r,r) = a(r)T 2 + c(r). This is just a limit case 
of the main curve. If we consider T + /3(ip) = (i.e. 
a = 0, c infinite) so that i/j = bfa, we see that the solution 
of Eq. (|4T|) is not unique for (3 = j3 mn since we can 
always superpose to Eq. (|C1Q[) an eigenfunction i/j mn of 
the Laplacian. Thus 



02 ->" 02 + X^n 



(F8) 



The condition T+j3(\jj) = 0, reducing to T/b+/3x(i>mn) — 
0, determines the value of x- Then, using Eq. (|67|) . 
we find that this solution exists for a unique value of 
the energy £? mn (r,r) given by Eq. (|F7[) . This is to be 
expected since c — » +oo when (3 — >• /3 mn . In conclusion, 
the case /3 — >> /3 mn or /3 = /3 mn with m and n odd is just 
a limit case of the main curve. 



3. Non-symmetric linear topography and r / T* 

The first case that will be treated is the one for which 
r ^ T* and T + p(i/>) ^ (i.e. a ^ 0, c finite). These 
solutions exist for any value of /3. The curve j3{E) relating 
their temperature to their energy is given by Eqs. (jUJ), 
(|39|) and (|45|) . This forms the main curve. Their entropy 
is given by Eq. (|49|) . For f3 = f3* we have c = leading 
to 1/E = 0. In that case, S/E = 0. 

Another possibility is T ^ Y* and T + /3(^) = (i.e. 
a = 0, c — » oo). In that case -0 = 602 which implies 
T + /3b (fa) = 0. Therefore, this solution exists for a very 
particular value of the inverse temperature that will be 
denoted /?o- The corresponding energy is given by Eq. 
([52]) where /? is replaced by fa. It turns out that this is 
just a particular point of the main branch corresponding 
to c — >• +oo so that it does not bring any new solution. 

In the case where /3 is equal to an eigenvalue f3 mn of 
the Laplacian, four cases must be distinguished. 

• Case 1: (3 = /3 mn with (Hi/j mn ) ^ and (Vw) 7^ 
(m odd and n odd) like /?n. In that case, Eqs. ([40]) 
and ([IT]) do not have any solution for f3 = f3 rnn and the 
functions fa and fa behave like 



, (Vw) , 

h ~ o _ o V>mn, 

P Pmn 



t>2 ~ ~- a Vmn, 



P-Pn 



(F9) 



(F10) 



when f3 — )• (3 mn . We are in the case T + (3(i/)) ^ (i.e. 
a/0, c finite). Substituting the asymptotic expansions 



([F9]) and (IfTO]) in Eqs. (|39j) and (|45j), we find that 

(Vw) 1 



(Hl/J mn ) Pmn(H^mn) 

I r 



and 



\b(Hl/) mn ) (l/jmn) 



(-. 



(p-p mn ), (Fll) 



/3 TOn \b(Hlp mn ) (ipmn) 



A, 



(F12) 



Therefore c and are finite as /3 — >> /3 mn . Substituting 
these results in Eq. ([48]) . we obtain a finite energy 



2£ 



1 



1 



(ipmn) 2 /3 n 



(H^ mn ) 
(i>mn) 



(F13) 



This solution is a particular case of the main curve. We 
note that the energy (]F13|) is of the form £? mn (r,r) = 
a(r)T 2 + b(r)T + c(r) with a > so that it forms a 
parabola. The minimum value of this parabola is E% 



reached for T = r™ n (r) with 



T--(r) (Ht/> mn ) 



J mn 



(Vw) 



(F14) 



• Case 2: fi = j3 mn with (H^ mn ) ^ and (ip mn ) = 
(m odd and n even) like /3i2. In that case, Eq. ([4T]) does 
not have any solution for f3 = /3 mn and 02 diverges like 



, (Hlpmn) 



(F15) 



when j3 — )• f3 rnn . We note that (02) = 0. Therefore, if 
T 7^ 0, we are in the case T-\-/3(ip) ^ (i.e. ft^O). Then, 
c is finite and (0|) -^ +oo so that 1/E — >> according to 
Eq. (HE]). If T = 0, we are in the case T + /3(ip) = (i.e. 
a = 0). Then, l/E->0 according to Eq. (|5^|). This is a 
limit case of the main curve. 

• Case 3: /3 = /3 mn with (Hi/j mn ) = and (Vw) = 
(m even and n arbitrary) like /?2i- The solutions of 
Eqs. (PU|) and (jUJ) are not unique since we can always 
superpose to Eqs. (|C6p and (|C1Q|) an eigenfunction , mn 
of the Laplacian. Thus 



01 ^ 01 + X^n 



02 ^ 02 + X^n 



(F16) 



(F17) 



The amplitude x ls determined by the energy constraint 
dm , replacing f3 by f3 mn and introducing Eqs. (|F16p and 
()F17p . This determines two types of solutions Xi(-B) and 
X2(E) for each value of the energy £^ (note that these 
two distinct solutions have the same value of c(E), hence 
the same value of a(E), given by Eq. (|45])). The en- 
semble of these mixed solutions form a plateau j3 = /3 mn 
starting from 1/E = (x — > +oo) and reaching the main 
curve at 1/E = 1/E mn (r,r) (x = 0). For 1/E = 0, 
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we find S/E — ^ /3 mn . On the other hand, E mn (r^r) is 
given by Eqs. (|35jL ([39]) and (|^5j) replacing /? by /3 mn and 
taking x — 0. It corresponds to a parabola of the form 
£mn(r,r) = a(r)r 2 + b(r)T + c(r) with a > 0. We shall 
denote by £"^ n (r) the minimum value of this parabola 
reached for T = r™^ n (r). The expressions of these quan- 
tities can easily be obtained from Eqs. ([48]) , ([39]) and ([45]) 
but they are not particularly simple so we shall not give 
them explicitly. Note that in a square domain (r = 1), 
since j3 mn = j3 nm by symmetry, we find that 1/E = for 
m even and n odd according to Case 2. In that case, the 
plateau reduces to a point. 

• Case 4: /3 = /3 mn with (H^ mn ) = and (Vw) 7^ 0. 
This is not possible for a topography of the form h = 
b(y - 2/0 ) with y 7^ 0. 



4. Non-symmetric linear topography and r = T* 

It will be first assumed that r = T* and T + /3(ip) ^ 
(i.e. a 7^ 0, c finite). These solutions exist for any 

value of /?. For /? ^ /?* , the curve /3(2£) relating their 
temperature to their energy is given by Eqs. ([48]) , ([39]) 
and ([^5]) . This forms the main curve. Their entropy is 

given by Eq. ([49]) . For f3 = /3i , the situation is different 
because Eq. ([^5]) takes in indeterminate form. In that 
case, c is determined by the energy constraint ([48]) with 

f3 replaced by /3* . This yields two solutions c\(E) and 
02(E) (and, correspondingly, two values of the chemical 
potential a±(E) and 012(E)) for each value of the energy. 
The ensemble of these solutions form a plateau /3 = /?* 
starting from 1/E = (c = 0, a = 00) and reaching 

the main curve at 1/E = 1/El (r) (c = 00, a = 0). 
When 1/E = 0, we find S/E = . On the other hand, 
^ fc) (r) is given by Eq. fl52|) with /? replaced by 0. 

A second possibility is Y = T* and T + f3(ip) = (i.e. 
a = 0, c — >- 00). In that case ^ = 6^2, which implies 
T* + /3b (fa) = 0. These solutions exist for j3 such that 
13 (fa) = /3*((fa)*)- Except in an antisymmetric domain 
where they exist for any f3 (see Sec. IIVA|) , in a non- 
symmetric domain they exist only for j3 = f3* and for 
the corresponding value of the energy E*(r). This is just 
a particular point of the main curve corresponding to 
c — » +00 so this case does not bring any new solution. 

The study of the eigenvalues f3 = f3 mn leads to the 
same results as in the previous section. 



Appendix G: Interest of the relaxation equations 

Relaxation equations for two-dimensional flows have 
been introduced by Robert & Sommeria [38j in the con- 
text of the MRS statistical theory and extended by Cha- 
vanis [26] to more general situations. These relaxation 
equations are relatively ad hoc since they are solely based 



on general "principles" similar to the first and second 
principles of thermodynamics: conservation of energy 
(and possibly other constraints) and increase of an "en- 
tropic" functional. Contrary to the Boltzmann equation, 
the relaxation equations are not derived from first prin- 
ciples but just assumed[45]. This is the reason why they 
exist in different forms. There is no general consensus 
about the importance and interest of these relaxation 
equations. Some researchers argue that they are useless 
and that they have no physical relevance. Other workers, 
on the other hand, consider that they are interesting for 
the following reasons: 

(i) They can be used as numerical algorithms to solve 
maximization problems of the form ([6]) with a general 
"entropic functional" of the form ([7]). This maximization 
problem can have various interpretations (as explained 
in Section Hi Bl and in [26]) related to thermodynamical 
and/or nonlinear dynamical stability of 2D Euler flows. 
Although the case of a quadratic functional ([24]) can be 
treated analytically (as done in our paper), the general 
maximization problem ©-([7]) is more complex and the 
relaxation equations can be a way (among others) to solve 
it. Indeed, it is generally difficult to directly solve the 
Euler-Lagrange equations associated with (j6])-([7j) and be 
sure that the critical points are entropy maxima. By con- 
trast, the relaxation equations guarantee (by construc- 
tion) that the relaxed state is an entropy maximum with 
appropriate constraints. Furthermore, when several en- 
tropy maxima (metastable states) exist for the same val- 
ues of the constraints, the relaxation equations can tell 
which (local) maximum will be reached by a given initial 
condition, thus unveiling its basin of attraction. Finally, 
by introducing in the relaxation equations a variable dif- 
fusion coefficient (related to the local fluctuations of vor- 
ticity), we can take into account the important effects of 
incomplete relaxation and lack of ergodicity [40|, ED] . This 
is not possible if we consider just equilibrium statistical 
mechanics. 

(ii) They provide non trivial dynamical systems, con- 
sistent with the equilibrium states, whose mathematical 
study is interesting in its own right. For example, we 
have observed that the system can remain blocked for 
a long time in an unstable state before finally reaching 
the maximum entropy state. Although this effect might 
be an artefact of the relaxation equations, the same phe- 
nomenon could also occur in more realistic parametriza- 
tions. Therefore, the study of the relaxation equations 
can suggest interesting ideas and developments. 

(iii) The discussion of Appendix [B] shows that the re- 
laxation equations are remarkably consistent with the 
parametrization of Holloway [36| based on relatively dif- 
ferent arguments. This is a hint that our parametriza- 
tion may contain ingredients that capture some features 
of real oceanic circulation. For example, Holloway has 
shown that the Neptune effect plays some role in the 
dynamics of the oceans. Therefore, making contact be- 
tween the two approaches is of interest. Furthermore, it 
is possible to include other geophysical effects in the re- 
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laxation equations such as wind forcing like in Kazantsev 



et al hfll. 
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